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THE RESPONSE OF AN AIRPLANE TO RANDOM ATMOSPHERIC DISTURBANCES 1 

By Franklin W. Diedeeioh 


SUMMARY 

The statistical approach to the gust-load problem which con- 
sists in considering flight through turbulent air to be a stationary 
random process is extended by including the effect of lateral 
variation of the instantaneous gust intensity on the aerodynamic 
forces. The forces obtained in this manner are used in dynamic 
analyses of rigid and flexible airplanes free to move vertically, 
in pitch, , and in roll. The effect of the interaction of longi- 
tudinal, normal, and lateral gusts on the wing stresses is also 
considered. 

The method of analyzing the rigid-body motions is similar 
to that used for analyses of the dynamic stability of airplanes, 
in that the equations of motion are referred to stability axes 
and expressed in terms of conventional stability derivatives. 
The method of analyzing the dynamic effects of structural 
flexibility consists in an extension of a numerical-integration 
approach to the static aeroelastic problem and is in a form 
which offers the possibility of calculating divergence and flutter 
speeds with relatively little additional effort. 

The mean-square values, correlation functions, and power 
spectra of some of the aerodynamic forces required in this type 
of analysis are calculated for one special correlation function 
of the atmospheric turbulence. It is shown, for instance, that 
if the span is relatively large compared with the integral scale 
of turbulence, the mean-square lift and root bending moment 
directly due to the gust are substantially reduced when the 
differences in instantaneous intensity of the turbulence along 
the span are taken into account. However, if the motions of 
the airplane are taken into account the mean-square root bend- 
ing moment may be increased as a result of these differences. 
Also, the meansquare pitching moment is shown to be sub- 
stantially increased if the tail length is relatively large compared 
with the scale of turbulence. Finally, the wing stresses due to 
longitudinal, normal, and lateral gusts are shown to be statis- 
tically independent under certain conditions. 

INTRODUCTION 

Tlio local velocity fluctuations acting on an airplane flying 
through atmospheric turbulence are functions of time defined 
only in a statistical sense and, hence, constitute a stochastic 
or random process. Consequently, the responses of the 
airplane, whether they are motions (linear or angular dis- 
placement, velocities, or accelerations), forces (lift, pitching 
moment, bending moment, and so on), stresses, or any other 


phenomena determined By the turbulence, can also he known 
as functions of time in only a statistical sense. 

This report is concerned with the statistical characteristics 
of those responses which have a bearing on the loads and 
stresses experienced by the airplane; although other prob- 
lems such as those relating to passenger comfort or to the 
stability of the airplane as a gun platform can be treated in 
the same manner, they will not be considered here. 

The first approaches to the gust-load problem which use the 
statistical techniques developed for stationary random proc- 
esses appear to be those of references 1, 2, and 3. (An earlier 
investigation concerned with the motions of an airplane in 
turbulent air is reported in ref. 4.) The fundamentals of 
these approaches are discussed in some detail in reference 1, 
and mention is made therein of investigations in other fields 
of engineering and physics that have dealt with the problem 
of deducing the statistical characteristics of the output or 
response of a dynamic system from those of its input. The 
mean-square normal acceleration of a rigid airplane free to 
move in one degree of freedom, namely, vertical motion, is 
calculated in reference 2. In addition to being rigid, the 
airplane is implicitly assumed to be small enough for all its 
components to experience the same gust velocity at any 
instant of time. This means that the span of the airplane 
must be small compared with the integral scale of atmos- 
pheric turbulence, which on the basis of the available 
knowledge concerning the properties of the atmosphere 
(ref. 4, for instance) appears to be in the order of several 
hundred to 1,000 or 2,000 feet; that is, the span of the airplane 
must be less than about 100 feet. 

The purpose of the present report is to extend this approach 
to large flexible airplanes free to move in all directions. As 
used herein, the terms “small” and “large” airplane refer 
to airplanes which are very small and not very small, re- 
spectively, compared with the integral scale of turbulence; 
thus, an airplane flying in a wide variety of atmospheric 
conditions may be “small” under certain conditions and 
“large” under others. Similarly, the terms “rigid” airplane 
and “flexible” airplane are used to designate airplanes flying, 
respectively, at speeds far below those at which dynamic 
and aeroelastic effects become important, and at speeds at 
which these effects have to be taken into account; the same 
airplane can thus be “rigid” under some conditions and 
| “flexible” under others. 


» Bupcrecdes NACA Technical Note 3910 by Franklin W. DIederich, 1937. This report represents, except lor some minor changes, a thesis submitted In May 19M In partial fulfillment 
of the requirements for the decree of Doctor of Philosophy at the California Institute of Technology, Pasadena, Call! 
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Several fundamental assumptions are inherent in the 
analysis contained in this report. In the first place, all 
atmospheric disturbances, motions, and structural deforma- 
tions are assumed to be small enough to produce forces that 
are linear and, hence, superposable. Also, the turbulent 
“input” to the airplane is assumed to be stationary in a 
statistical sense; that is, the turbulence in the plane of the 
flight path is homogeneous. For the large airplane, the 
additional assumption is made that the turbulence is usd- 
symmetric with respect to vertical axes, a condition less 
severe than complete isotropy. The statistical character- 
istics of the turbulence are thus assumed to be invariant 
under a translation of the space origin within the horizontal 
plane and under a rotation of the coordinates about the 
vertical axis. Finally, Taylor’s hypothesis to the effect 
that time displacements are equivalent to longitudinal space 
displacements is assumed to be valid. 

The aerodynamic forces directly due to atmospheric 
turbulence, which constitute the input forces for the dynamic 
system represented by the airplane, are calculated in the 
first part of this report for the large airplane, that is, for the 
case where the spanwise distribution of the intensity of 
turbulence has to be taken into account. (The effect of 
spanwise variation of gust intensity on the lift has been 
treated by a slightly different method in ref. 5.) 

The dynamics of the rigid airplane are considered in the 
second part. The dynamic system is now represented by a 
set of three simultaneous ordinary differential equations, 
rather than one as in reference 2; nonetheless, the problem 
of calculating the required transfer functions is still one of 
simple algebra. 

The next part is concerned with the small flexible airplane 
and thus has direct application to fighter-type airplanes and 
guided missiles operating at relatively high speeds, in addi- 
tion to serving as a preliminary to the last part of the report. 
The dynamic system is now represented by a partial differ- 
ential equation, and the calculation of the transfer functions 
requires the solution of ordinary differential equations. 
Once these functions are calculated, however, the statistical 
techniques are the same as before, as a result of the fact that 
the lateral variation in gust intensity is ignored. Either 
modal or numerical-integration approaches may be used to 
analyze the dynamics of a swept-wing airplane with arbitrary 
stiffness and mass distribution. Although modal approaches 
have usually been preferred in the past for similar problems, 
it was believed that, in view of the highly complex nature of 
modern aircraft structures and the advanced type of com- 
puting machinery required and generally available for their 
analysis, the numerical-integration approach would be 
preferable and it has, therefore, been used. 

The last part contains the analysis of the large flexible 
airplane. The statistical problem is now that of a system 
which is characterized by a partial differential equation 
with time and a space coordinate as independent variables 
and which is subjected to a random input that varies in time 
and space, so that more is required than the transfer func- 
tions from the gust intensity at one point on the wing to 
the stresses at another. The particular statistical problem 


presented by this case is considered in some detail, and the 
appropriate transfer functions are then obtained by using 
the numerical-integration approach presented in the preced- 
ing part to solve, in effect, the ordinary differential equations 
which describe the wing deformations at any given fre- 
quency. 

SYMBOLS 

A aspect ratio 

b span 

G(k) Theodorsen function 

C L . lift coefficient, LfqS 

G l lift-curve slope 

O t rolling-moment coefficient, L'/qSb 

Oi coefficient of damping in roll, defined as positive 

for positive damping 
G m<[ damping-in-pitch derivative 

C Ma static pitching derivative 

c chord, parallel to plane of symmetry 

c average chord, S/b 

c t section lift coefficient at station y, l/qc 

El bending stiffness 

e } distance from section aerodynamic center to shear 

center, fraction of chord 

es distance from shear center to section center of 

gravity, fraction of chord 

e 3 distance from shear center to the midchord point, 

fraction of chord 

e t distance from shear center to the %-chord point, 

fraction of chord 

Q(y, 17) dimensionless lift-influence function (Green’s func- 
tion for the spanwise lift distribution) 

GJ torsional stiffness 

g acceleration due to gravity 

H(o>) response to sinusoidal oscillation, Fourier trans- 

form of h(t) 

h(t) indicial-response function 

I x mass moment of inertia about X-axis 

I v mass moment of inertia about F-axis 

Jo, Jr Bessel functions of the first kind, order 0 and 1 

Kq,Ki modified Bessel functions of the second kind, order 
0 and 1 

Kiu integral of Kq 

k reduced frequency, a>c/2 U 

k' dimensionless frequency, uL*jU 

L lift 

L' rolling moment 

L* integral scale of turbulence 

l distributed lift per unit distance along the span 

M pitching moment 

M b bending moment 

M t twisting moment 

m mass (of airplane, unless designated otherwise by 

subscripts) 

m distributed mass per unit distance along the span 

m T distributed twisting moment (about axes perpen- 

dicular to the plane of symmetry) per unit 
distance along the span 
q dynamic pressure 
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r radius of gyration about center of gravity; longi- 

tudinal displacement corresponding to time 
displacement r 

5 wing area 

s scale parameter, 

t time 

U mean flying speed 

u longitudinal component of gust velocity 

u p horizontal component of disturbed motion 

v lateral component of gust velocity 

W weight of airplane 

w vertical component of gust velocity 

w p vertical component of disturbed motion 

x coordinate along mean flight path 

Ax. distance from intersection of elastic axis and root 

chord to airplane center of gravity 
x, tail length, distance from airplane center of 

gravity to aerodynamic center of tail 
x{ modified tail length, distance from intersection 

of elastic axis and root chord to aerodynamic 
center of tail 

y coordinate perpendicular to plane of symmetry 

z coordinate ih plane of symmetry perpendicular to 

mean flight path; vertical deflection 
a inclination of chord to X-axis 

/3 span ratio, b/L* 

T (y) autoconvolution function for y(y) 

y(y) dimensionless lift distribution, ccJcOl 

y' (y) dimensionless lift distribution in roll, cc,[cCi r 

y variable of integration corresponding to y 

6 angle of pitch 

A angle of sweepback 

p mass density of the air 

r time displacement, argument of time-correlation 

function 

<f> the power spectrum $> in the case of axisymmetry 

'T> two-dimensional power spectrum (double Fourier 

transform of ip) 

<p(k) Sears function (unsteady-lift function for gust 
penetration) 

i v one-dimensional or point power spectrum 

<p two-dimensional power spectrum (single Fourier 

transform of ip) 

P one-dimensional or point correlation function 

P two-dimensional correlation function 

co frequency of oscillation 

Subscripts: 

<?,e effective 

/ fuselage 

r wing root 

t tail 

u horizontal component of turbulence 

w wing; vertical component of turbulence 

Matrix notation: 

[ ] square or rectangular matrix 

l 1 diagonal matrix 

L J row matrix 

{ } column matrix 


Dots over symbols indicate derivatives with respect to 
time. 

AERODYNAMIC FORCES RESULTING DIRECTLY FROM 
ATMOSPHERIC TURBULENCE 

The motions of a rigid airplane depend on the overall 
forces and moments, whereas the stresses of a rigid airplane 
and the motions and stresses of a flexible airplane depend on 
the distribution of these forces, as well. This part of the 
report is concerned with the calculation of the integrated 
and distributed forces and moments directly due to atmos- 
pheric turbulence when the spanwise variation of gust 
intensity has to be taken into account. (The forces and 
moments caused by the motions which result from the forces 
treated in this part can be calculated by conventional 
methods and will not be considered here, although the com- 
bined forces will be considered in the following parts.) 
Thus, this part serves as a basis for all the material presented 
in the later parts pertaining to the large airplane, and, hence, 
the fundamental notions required for an analysis of the large 
airplane are introduced here and discussed in some detail. 

The basic approach is as follows: First, the instantaneous 
value of the quantity of interest, such as the lift, is expressed 
in terms of the instantaneous gust intensity at a point and a 
suitable influence function. For the lift this influence func- 
tion can and will be identified with a certain lift distribution 
on the given wing in reverse flow; the same procedure may 
and, in the case of the rolling moment, will be followed for 
other integrated forces. On the other hand, for the local lift 
the influence function is the Green’s function for the three- 
dimensional unsteady-lift problem and cannot be identified 
with an easily calculated lift distribution on the wing in 
reverse flow. Inasmuch as no knowledge concerning this 
function appears to be available, a method of calculating an 
approximate Green’s function for this problem is outlined 
herein. The required influence functions for integrated 
effects can be synthesized from this function, and if the 
associated lift distribution in reverse flow cannot be calculated 
conveniently this approach may be preferable. This tech- 
nique is illustrated here by means of the bending moment. 

The next step consists in using the expression for the 
instantaneous value of the given quantity to calculate a 
correlation function for this quantity in terms of a correla- 
tion function of the normal component of the atmospheric 
turbulence. The power spectrum for the given quantity 
can then be obtained by taking the Fourier transform of 
its correlation function. This power spectrum is considered 
herein to be the desired end result, because the mean- 
square values of the quantity and its derivatives can be 
obtained from it, and other statistical parameters of interest 
can be obtained from these mean-square values. Several alter- 
native approaches for calculating the aforementioned spec- 
trum, either from the correlation function or directly from 
the spectrum of atmospheric turbulence, are given in con- 
nection with the lift and are directly applicable to other 
quantities as well. 

In this part of the report the assumption is made that the 
influence functions of concern can be written as products 
of a function of time alone and a function of distance along 



146 


REPORT 1345 — NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


the span alone, and advantage is taken of this simplification 
in calculating the desired spectra. This restriction is 
removed in the last part of the report, where the influence 
functions considered cannot be separated into space- 
dependent and time-dependent constituents, and so the 
approaches outlined there are generalizations of those 
presented in this part; they may also be used for the calcula- 
tion of the spectra of the quantities considered in this part 
if the assumption concerning the influence functions is 
not valid. 

DEFINITIONS OF STATISTICAL PARAMETERS 

As pointed out in the introduction, the intensity of the 
vertical component of turbulence w(t) is a random process, 
so that the resulting forces can also be known only in a 
statistical sense. The purpose of this part is to calculate 
certain statistical properties of these forces, namely, their 
mean-square values, their correlation functions, and their 
power spectra. The fundamental principles involved in 
statistical analyses of the type considered herein are ex- 
pounded in some detail, and citations of the literature on 
the subject are given in references 1 and 2. These funda- 
mentals will therefore not be repeated here. However, 
both for the sake of ready reference and inasmuch as the 
statistical terms are not always defined in the same manner, 
the forms that are used herein are indicated in the succeeding 
paragraphs. 

The time average of a time-dependent quantity is design- 
ated by a bar placed over the symbol, 2 and is defined as 
follows: 

i- lin x m f Z(t)dt 

T->™ 4J- J-T 

The assumption will always be made that this limit exists 
and is invariant under a translation of the origin of time. 
This assumption implies that the processes considered here 
are stationary in a statistical sense. 

The mean of a random process /(f) is defined as its time 
average, and is always assumed to be zero. In dealing 
with processes with nonzero mean this analysis is thus 
pertinent only to the process which consists of the difference 
between the original process and its mean value. Similarly, 
the mean-square value of a random process /(f) is defined 
as the time average of the square of the process, so that 

y*(f )dt 

T-*<» A-L J-T 

The time-correlation function of /(f) is defined as 

$Ar) 

so that 

T'=M0) 

and the power spectrum of /(f) is defined as the Fourier 

3 When no possibility of confusion exists, a bar la also used to designate a space average, 
as in the case of c and jr. Also, for the components of turbulence the correlation functions 
depend primarily on space displacements and can be defined by space averages. For the 
s ake of consistency, however, they are considered to be defined by time averages. 


transform of the correlation function : 



The second form is the more convenient one when repre- 
sents more nearly a space correlation than a time correlation, 
so that it depends directly on a space displacement rezUr 
rather than on a time displacement r. For the purposes of 
the analysis presented herein, the gust or input correlation 
functions have this property, so that, for instance, may 
be defined as 


Wd=w( t)w(t+^ 

By virtue of the reciprocal properties of Fourier trans- 
forms and the symmetry of //r), which is a direct conse- 
quence of the assumed stationarity of /(f), J 3 may be ex- 
pressed in terms of the spectrum of / as 

**’)*' (1) 

If /(f) is the input of a linear system, the power spectrum 

of the response x(t) of the system is related to the power 
spectrum of /(f) by 

ft(«)Hff(«)|V/(») (2) 

where H(o>) is the transfer function of the system, that is, 

the complex amplitude of the response of the system to unit 
sinusoidal input. Hence, H{co) is also the Fourier transform 
of the indicial response h(t) of the system, which is defined 
herein as the response of the system to a unit impulsive 
input: 

H(«)= r e~ ia, h(t)dt (3) 

J — ao 

where the lower limit could be taken as zero, since /(.(f) is 
zero for f<0. Conversely, h(t) can be obtained from H(w) 
by means of the inverse of equation (3) : 

h(t) e ta ‘H(u>) dco (3a) 

The mean-square value of the response can then be ob- 
tained by integrating its spectrum. Similarly, the mean- 
square values of the nth time derivative of the response can 
be obtained from the (2n) th moment of the output spectrum. 
For instance, 

From the mean-square values of these derivatives other 
statistical quantities of interest, such as the expected number 
of peaks of the response per unit time, can then be calculated. 
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LIFT-INFLUENCE FUNCTIONS IN UNSTEADY FLOW 

At any time the lift on a wing which results directly from 
atmospheric disturbances can he expressed for an unswept 
wing as 

rw pin 

A(t)== dh h{ti,y)w(U{t—t^,y)dy (4) 

J-a J-W i 

where h(t,y)dy is a lift-influence function which represents 
the lift caused by an impulsive vertical gust of width dy 
which at time f=0 impinges on the wing at station y. 

The influence functions required in equation (4) are diffi- 
cult to calculate directly; methods for obtaining lift distri- 
butions on wings of finite span in unsteady flow usually 
require numerical solutions which do not lend themselves 
readily to the analysis of angle-of-attack distributions repre- 
sented by delta functions. However, by virtue of the reci- 
procity theorems of linearized lifting-surface theory (ref. 6, 
for instance) the lift influence function for a twisted wing in 
indicial motion is equal to the lift distribution on that wing 
during indicial motion in the reverse direction with uniform 
unit angle of attack. The lift distribution in indicial motion 
with uniform angle of attack can be calculated relatively 
easily. 

For the few cases for which calculations have been made 
(namely, some unswop t wings), this lift distribution tends to 
be substantially invariant in time, except for overall magni- 
tude. For instance, the calculations of reference 7 indicate 
that the lift distribution of an oscillating rectangular or 
elliptic wing in incompressible flow is substantially inde- 
pendent of frequency, so that in indicial motion it is sub- 
stantially independent of time. This simplification may not 
be valid for swept wings. 

The lift influence function can then be written as 

h(t>y)=\^t) T'6/) (5) 

where y{y) defines the steady-state lift distribution for 
uniform unit angle of attack: 



and where h L (t) describes the variation of the overall magni- 
tude of the lift as a function of time after entry into a sharp- 
edge gust and may be written as 


h L (t) 


d 

U dt 



In turn, k 2 is the lift response to a unit sharp-edge gust 


normalized to a steady-state value of unity and is, as in the 
preceding equation, usually expressed in terms of semichord 

lengths traveled The time derivative is taken here be- 
cause the response wanted is the one to a unit impulsive 
gust rather than a sharp-edge gust. 

The Fourier transform H L (u) of this function h L (t) is 
proportional to a function <j>(k), which may be termed the 
generalized Sears function because for two-dimensional in- 
compressible flow it is the Sears function: 

(6) 

where 



The function H L (o>) represents the complex amplitude of the 
lift due to sinusoidal gusts of unit amplitude and is thus the 
transfer function from the gust to the lift. 

Actually, only the absolute square of 4>{k) will be required. 
For two-dimensional incompressible flow the following 
approximation is given in reference 1 : 

i«»r-rvh < 7) 

This expression has the advantage of simplicity, although it 
is somewhat in error compared with the absolute square of 
the Sears function at very low frequencies, a fact which could 
be remedied by using the approximation 

1^® l 2ta l+aia+2^). 

where a is about 15 for a good overall fit to the exact expres- 
sion. However, the behavior of the Sears function itself at 
very high frequencies is unrealistic, because its absolute 
square goes to zero as 1/k, whereas for any nonzero Mach 
number and any finite span the absolute square of the gen- 
eralized Sears function can be shown to tend to zero at least 
as rapidly as l/P. As a result of these discrepancies, the 
approximation given by equation (7) and, for the same reason, 
the absolute square of the exact Sears function, cannot be 
used to obtain moments of the lift spectrum, that is, values 
of the mean-square derivatives of the lift, although they may 
be adequate for calculating the mean-square value of the 
lift itself in many cases, particularly when absolute accuracy 
is not required. Whenever possible the values of |$(fc)|* 
given in reference 8 for the plan form and Mach number of 
interest should be used. 


MEAN-SQUARE LIFT AND ITS SPECTRAL RESOLUTION FOB THE UNSWEPT WING 


Basic equations. — The correlation function of the lift can be expressed, by virtue of equation (4), as 
/■ * c« Pin pin : 

«fc(V)= h(ti,ydh(t i ,ydw<JJ{t—t l ),y 1 )v)(!J(t-)rT—U),y 3 )dy l dyylt l dU (8) 

J -°= J -u 2 J -Hi 

where the averaged product on the right side represents a velocity correlation function. This function depends in general 
on hoth space and time displacements. However, if Taylor’s hypothesis is made, the time displacements are equivalent to 
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longitudinal space displacements. The velocity correlation functions are then functions only of longitudinal and lateral 
space displacements. Thus, for homogeneous turbulence, 

w(x+Ut,y)w(x+^-\-U(t+T),y+y)=^ w (^+UT,y) (9) 

/ 

In addition to Taylor’s hypothesis and the assumption of homogeneity, the turbulence is assumed to be axisymmetric 
with respect to vertical axes, so that is a function only of V? 2 + 7 ? 5 - This function is the ordinary space-correlation 

function $ a , so that 

$ a (U V+f, y) =^„(V (U r+£) 2 ++) 

and, hence, 


rt>n rm . 

h(tuy^h(ta,y^ a b/U 2 (T+ti—t 3 ) i +(y i —yi) i )dyidy : dt 1 dti 

-■» J—a J ~lf 3 J-b/2 


If the assumption imphcit in equation ( 6 ) is now made, 
the preceding equation can he written as _ 



h L (ti) hzfa) 'Pa^QJ (r — (— i i — fj)) dti dti 


( 10 ) 


where 

^ B (y t/ J T s + {rji—yi) 2 ) Y(yi) y(y*) dy I dy 2 
=^J' o T(tj) twUU'r'+r?) dr] (11) 


where, in turn, Y(y) is an autoconvolution of y(y) defined by 

b 

r(’))=|J_^ y(y) y(y+v) dy (12) 

The validity of equation ( 11 ) can he demonstrated by 
performing the integration in the y u y 3 plane as indicated in 
the following sketch: 



That is, integration is performed first over y u with a variable 
7 j =y 3 —y 1 held constant. For this integration is constant, 
so that only 7(1/1) and 7(2/1+17) are involved, and the result 
is a function of 17 which is one-half of the function r(i?) de- 
fined by equation ( 12 ). The second integration is then per- 
formed over 17, yielding equation (11) except for a factor of 


In this process, only the part of the square above the lino 
2 / 2 = 2 /i is covered. However, by a similar process, the part 
of the integral corresponding to the part of the plane below 
this line can be evaluated and shown to be equal to the first 
part, so that by defining T as in equation ( 12 ), both parts 
are taken into account simultaneously in equation ( 11 ). 

The quantity 


=jJ 0 ^»( 1 ?) dy 

may be considered an averaged mean-square vertical com- 
ponent of turbulence; ^„,((7 t) is then the corresponding cor- 
relation function, and the Fourier transform of the latter, 

^>)=iJV^ (Pr) f Kt (Ur) d(Ur) (13) 

is the corresponding power spectrum. 

Once <p w ,(oi ) has been obtained, the power spectrum of the 
lift can be obtained by taking the Fourier transforms of both 
sides of equation ( 10 ). The result is 

%.(w)=|fli(w)|Vu,(«) (14) 

where H L {u>) is the transfer function defined in the preceding 
section. The mean square of the lift can then be obtained 
byintegrating its spectrum, as indicated in equation ( 1 ). 

Equation (14) has the same form as the corresponding 
equation for the case where spanwise averaging of the effects 
of turbulence is not taken into account (see eq. ( 2 ) and ref. 
1 ), except that <p K (u>) is now replaced by <p Ut (o>). Thus, the 
spectrum of the averaged turbulence must approach that of 
the unaveraged turbulence when the span approaches zero, 
as may be seen to be the case from equation ( 11 ) and the 
definitions of 7 (y) and T(y). 

Two alternative approaches. — The defining relations for 
equations ( 11 ) and (13), do not necessarily represent 
the best method of calculating it in any given case. A sligh tly 
different expression appears to be more convenient in general. 
It consists in substituting equation ( 11 ) into equation (13) 
and inverting the order of integration to yield 

<Pw t (a)=2; r(y)p w (co,y) dy (15) 



THE RESPONSE 0E AN AIRPLANE TO RANDOM ATMOSPHERIC DISTURBANCES 


149 


where 

?»(<», v)=^j j e t?)d(U r) . (16) 

bo that 

<pw(w,0 )=¥>»(w) 

A third approach, which has certain advantages over the others, is similar to the one that has been used in reference 5. 
In this approach the assumption of axis ymm etry is not made initially, and use is made of the spectrum 

5 w (Xi,X 3 )=^ f“ dz dr, (17a) 

7T J — CO J — CD 


The correlation fimction ^ tt> (Z7V) can then be written as 


i nt/i rb/i i ra> r a _ 

J J y(y l )y(y 2 ) dy, dy a £ J_ = ^(X,,^) d\, d \ 2 


Substitution of tliis expression into equation (13) and then 
interchanging the order of integration yields : 

^ ( “ )= 2(y/l 5^>x 2 )|r(x 2 )| 2 <& 3 

where 

~ 1 r bn 

r (X)=r y(y) e-^ dy 

Vj-blZ 

and is, as a result of the symmetry of y(y), real and sym- 
metric in X. 

Now, if the turbulence is axisymmetric, depends only 
on Xea^Xi s -t-X 3 a , that is, 


? K (X,,X 3 )=$„(VX?+V) (17b) 


so that the expression for <p We { «) becomes 

J 0 ” $, {V(5) i+X1 )^ (X)<ZX (18) 

where $„(X) is related to ^„(r) by 

$o,(X) =- f r Jo(Xr) #„(r) dr 
TTjo 

Although the spectra used in the last two approaches have 
been defi n ed by expressions involving the point correlation 
function they can be expressed equally well in terms 

of the point spectrum <p w (u>) , so that if, say, an experi- 
mentally obtained point spectrum is to be used, it need not 
be transformed into a correlation function before it can be 
used in these calculations. The required relations are 

Vw(a,y)=<pw(w)~- SjjJ^Pafo) -j=== Ji tftt! 

(19) 

and 


where the notation f~ is used to specify that the finite part 
of the integral is to be taken, an operation which may be 
performed by integrating by parts and ignoring the infinite 
part, so that, in terms of a proper integral, 


*„(X) 


2TJ 2 fi" d<p w (id) dw 
TT Jxt7 d(0 


Also, the function $„(X) can be obtained from v a (oi, y) by 
means of either of the relations 


$»(X)=JyJ' e i * , p ttl (0,y)dy 

2 r°° - 

==^J VwQdJ,y)dy 


Results of calculations. — In order to illustrate the magni- 
tude of the effects under consideration, calculations have 
been made for a uniform loading y(y)=l and a point correla- 
tion function which has been used in references 1 and 2 and 
appears to fit experimental data (ref. 4) fairly well over a 
large portion of the significant frequency range, namely, 

so that 


i 


<pa(cc) 


u?L* l+3k' 2 
: idJ (!+£'*)* 


where L* is the integral scale of turbulence, which is here 
defined as 


£*== f fw(r) dr 

w*J o 

and where 



$»(X)=— 


2U 2 

7T 



CO 

(co s -X 5 Z7*)« 


du> 


( 20 ) 


This correlation function has the drawback that the mo- 
ments of the spectrum associated with it are infinite, so that 
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it implies a process with infinite mean-square derivatives, 
but it is quite useful if only the mean-square value of the 
process itself is of interest. 

For uniform loading, 


and 


r(,)=2(i— |) 


sm ■ 


\6 


r(x)= 


\b 

2 


and for the given correlation function, 

\j* (T+IV^ 1 (x^ /T + Fl )- 

(z*) lTF^ 0 (^ /I + Fi )] 


where Kq and Ki are modified Bessel functions of the second 
kind, and 


$«(X) 


ZvfL** fc'* 

(1+F 1 ) 54 


The mean-square average gust intensity for this case is 
given by 


w,*— w* 


1— e~* 


where fi^b/L*, and is shown in figure 1 (a). The correlation 
function is given by 


[*&(*; sinh- 1 ^)-^^; sinh^-h 



where a^r/L*, and where Ko{<r;d) and Ki(<r;d ) are incomplete 
modified Bessel functions of the second kind defined by 

i£„(a",d) = j cosh nd e~ a 00311 0 dd 

The spectrum <p„ e (&)) is given by 

¥>*»=! ^ [KM T+F 5 )- 

+ (1-3*' 1 ) [2- 

2ejiw*K 1 (jhtT+k*) -fpa+k'^Koihfiw 2 )]} 

where Kio{x) is the integral of Kq(x) : 

Zio(x)=£ Koix^dxi 

The functions \f/ Wt and <p„ e (normalised with the averaged 
mean-square turbulent velocity) are shown in figures 2 and 




(b) Swept wings, — =0.6. 

L* cos A 

Figure 1. — The mean-square averaged vertical component of 
turbulence. 

3, respectively. The effect of the span ratio on the normalized 
correlation function of the averaged turbulence may be seen 
to be relatively small, so that the effect on tho unnormalized 
correlation function is primarily the decrease in overall 
level given hy the ratio of w/ to vf. A similar statement 
may be made for the power spectrum. If the power spectrum 
were not normalized the averaging effect of the span would 
tend to reduce the intensity of the spectrum at all fre- 
quencies, but the high frequencies would be attenuated 
much more than the low ones, as might be expected. In 
fact, although the unaveraged spectrum decreases as o>~ 2 at 
high frequencies, the averaged spectrum decreases as w 3 . 

The asymptotic values for b/L*->co are shown in figures 
2 and 3 in order to indicate the nature of tho functions 
considered here when the scale of turbulence is small com- 
pared with the span, as may be the case for a wind-tunnel 
model responding to natural or artificial tunnel turbulenco, 
or for a buffeting wing or tail surf ace, although tliis condition 
is not of practical concern for the gust-load problem. 

The power spectrum of the lift is equal to the product of 
the power spectrum <pa,(u>) and the absolute square of tho 
transfer function S L (cd), as indicated in equation (14). 
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Figure 2. — The correlation function of the averaged vertical 
component of turbulence. 


integrating the spectrum. Thus, for instance, the approxi- 
mate expression for [<£(&) | J given by equation (7) and the 
spectrum <p a (o>) used for the preceding calculation yields 
.the mean-square lift: 






Figure 3. — The power spectrum of the averaged vertical component of 

turbulence. 

Inasmuch as this lift is not an end in itself but only one 
of the parameters that enter into the calculations of the 
motion of the airplane, its mean-square intensity is of little 
practical significance; its spectrum is the quantity needed 
in further calculations. However, if the mean-square 
intensity is wanted for any reason it can he obtained by 


and, in view of the observation that much of the turbulent 
energy is contained in a region for which the span has a very 
small effect on the (normalized) spectrum, this equation 
should serve to furnish an approximation to the mean- 
square lift for nonvanishing span, provided w t * is used 
instead of w 2 . However, the mean-square values of the 
derivatives of the lift cannot be calculated in this simple 
manner, because the deviation of the normalized averaged 
spectrum from the point spectrum at high frequencies cannot 
be ignored in calculating the moments of the lift spectrum. 

MEAN-SQUARE LIFT OF THE SWEPT WING 

For the yawed or sideslipping unswept wing, equations 
(10) and (14) for the lift-correlation function and spectrum 
are still valid if an appropriate lift-influence function is used, 
and if the correlation function for the averaged turbulence 
is defined by 

’Aw, (U r) — A J o r (y) 'l'wU(Ur+ 17 tan A)M-7?*) dy 

where T(y) now pertains to a lift-distribution function 
y(y) which is appropriate for yawed motion and is defined for 

— — cos cos A. The mean-square averaged intensity 

of the vertical component of turbulence is then 


Wt *= ■■■■— p A T( v ) fJ— j) dy 
b cos A Jo ^cos A/ 

=£ J* o Tfij'cos A) \A k (ijO dy' 


where y' = — —• Thus, this mean-square intensity is un- 
cos A 

affected by the yawing process, except for the slight change 
which results from the change in T(y), although the spectral 
resolution of the averaged intensity changes in the process. 


For the swept wing both 2/2 — Vi and |yj| — |i/i| occur in the integral, so that the reduction of the double integral for 
'I'u^Ut) to a single integral (see eq. (11)) cannot be affected so simply. The double integral for the swept wing is 

AUt)=± 7 ( 1 / 1 ) y{y») tu UlUr- Klifel-ty.1) tan Al^C^-y,) 2 ) d Vl dy 2 

0 J -mJ — */a 


G2C007 — 01 


•11 


(21) 
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From this integral, by using rectangular lift distributions 
and the aforementioned point correlation function, the mean 
square of the averaged turbulence as well as the corresponding 
correlation function and power spectrum have been calculated 
by numerical integration for various sweep angles A, the 
0 

ratio r being maintained at 0.5. (The decision to hold 

cos A & v 


P 


■ rather than /3 constant was reached as a result of the 


cos A 

foregoing analysis of the yawed unswept wing, which indi- 
cated that the effects of sweep should be minimized in this 
manner.) The results for w , s are shown in figure 1, and 
the effect of sweep on w S s is seen to be small for this com- 
parison. The calculated correlation functions and spectra 
(normalized with respect to w,*) are not shown because they 
agreed with those for A=0 within less than 1 percent for 
most values of Ut and k', respectively. 


MEAN-SQUARE ROLLING MOMENT 

In the preceding sections the averaging effect of the span 
has been shown to consist, essentially, in reducing the 
effective intensity of the turbulence sensed by the wing; thus, 
it only modifies the forces present on a -wing of small span. 
If the analysis is extended to the rolling moment, however, a 
new phenomenon appears. When a wing is so small relative 
to the scale of turbulence that at any instant all of its points 
experience the same turbulent velocity, the wing experiences 
no rolling moment as the result of the direct action of 
turbulence (although it may experience a rolling moment 
indirectly as a result of the rolling and yawing motion caused 
by the lateral component of the turbulence) . On the other 
hand, on a large wing the different intensities of the turbu- 
lence at different points on the span give rise directly to a 
net rolling moment, which then results in rolling motion. 
In this section the mean-square value of this moment is 
calculated. 

At any instant t the rolling moment IJ (f) can be written 
in the same form as the lift L(t) in equation (4) ; however, 
according to the previously mentioned reciprocity theorem, 
the lift-influence function h(t,y) is now the lift distribution 
for an indicial roll with unit helix angle at the wing tip. If 
the assumption of invariance of this distribution with time 
is made, as for the symmetric case (see eq. (5)), then the 
required lift-influence function can be written as 




Ci(t)gSb 

where h'(t)= ^ ^jj — , and where the steady-state lift 

distribution y'(y) 1 ==^- now pertains to a unit linear anti- 
symmetric angle of attack. 

The correlation function for the moment can then be 
written as 

J hf (tdh' t^))dtidti ( 22 ) 


where 

'Pw c '(Ut)=^ f f Y (y,) y' (yMw tiU*T*+ (j/s— 2 /i)*) dy x dy, 

o J -bja J -b/2 

(V Y)dv (23) 

where, in turn, in analogy with equation (12), 
o riw)- n 

r'0»)=jj t '(y)y'(y+v)dy 

Hence, the second and third approaches indicated in the 
section concerned with the mean-square lift (see eqs. (15) 
and (16)) can be used to obtain 

1 P 

¥>«,,' («)=£ I T'(v)<Pw(oo,v)dy (24) 


and 

(V(&T +XJ )l f/(x) l yx 

where 

_ 1 rbl 2 

r '(>0=£ J ^y' (y)e~ i)Jl dy 

and 


gff 1 ) \<t>(k)\ 2 <p w ;(to) (26) 

so that the mean-square rolling moment can be obtained by 
integrating this spectrum. 

A qualitative indication of the effect of span on the mean- 
square rolling moment sensed by an airplane may be obtained 
from the quantity ’/'«■/ (0), which represents the integral of 

the spectrum <pv,'. For a linear loading, and 


r' 





r / (X)=6i 


\b \b 
— cos -sin 


U 

2 


Hence, using the aforementioned expression for the point 
correlation function yields 


T "1 Q 

«V(«) { [-32+20'*+ (l 6/3 '‘+^)W) + 

(32p'+Gp' 3 )K 1 (P')]+k' i [32-Qp' i -m' t K 0 (p')~ 

(32f}'+2fJ'*)K 1 (P')+p' 3 Ki 0 (f}')] } 


where 


p^/sVI+F 5 
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and 

^'(°)= i |?[f-6+^^(6+6^+3^+|/3 3 )] 


from elementary definitions and from the aforementioned 
reciprocity theorem: 



This spectrum and its integral approach zero linearly as 
tends to zero. 



yn(y) 0i(;y) dy 


GENERALIZED AERODYNAMIC INFLUENCE FUNCTIONS IN 
UNSTEADY FLOW 


The aerodynamic influence functions used in the preceding 
sections define the contribution of a given station of a wing 
to the total lift and rolling moment. In the analysis of a 
flexible wing, and even in the calculation of certain properties 
of a rigid wing, generalized aerodynamic influence functions 
are required, which define the contribution of one station on 
the wing to the lift at another station and thus represent a 
Green's function for the unsteady spanwise lift distribution. 
No work appears to have been done on such functions. For 
steady flow, apart from some calculations for supersonic 
speeds which are based on the subdivision of a given wing 
into a number of squares, the only available results appear 
to bo those given in references 9 and 10. 

The analysis in this section is based on reference 8 and 
consists in a generalization of the method presented therein 
to unsteady flow. This method constitutes an attempt to 
predict the lift distribution for any given twist on the basis 
of knowledge concerning a few definite angle-of-attack dis- 
tributions, and may therefore be termed a function-interpo- 
lation method. For the present purpose, the presumably 
known lift distributions are the ones for uniform angle of 
attack in direct and reverse flow in a dimensionless form, 
namely, cc,/cC L ; they will be referred to as y D (y) and y s (y), 
respectively. (The function y(y) used previously is the one 
now designated by y Also required is the coefficient 
of damping in roll (7j p . (The lift-curve slope and the coeffi- 
cient of damping in roll are the same in direct and reverse 
flow by virtue of the reciprocity theorem, so that no distinc- 
tion will be made.) 

The approach of reference 9 then yields the following 
approximate expression for the lift distribution due to any 
angle-of-attack distribution: 


^=C La {Sy D (i/) + K\a(y)—a\ y R (y)} 

where 


l rb / s 

a= bJ V2 7fl ^ dy 


K=-°' v 


b 3 


r-bn 

“ 4 Jo 7 *^ yid y 


(26) 


(27) 


Values of K may be obtained from the information given in 
references 9 and 10. As the aspect ratio tends to zero, K. 
approaches ]i, whereas for aspect ratios approaching in- 
finity, K tends to 1. The following relations can be obtained 


l rbn 1 rbn 

fcj-ws dy= bJ m 7r(v ‘ I dy=1 




With the aid of these relations the lift distributions given 
by equation (26) may readily be seen to have the correct 
lift and rolling moment for all angle-of-attack distributions, 
and to reduce to the exact lift distributions for angle-of- 
attack distributions which vary linearly along the span. A 
lift distribution which possesses these properties could readily 
be obtained by approximating any given angle-of-attack 
distribution by a suitably chosen linear one. However, this 
procedure would yield poorer approximations than the ones 
furnished by equation (26); for a parabolic angle-of-attack 
distribution on a wing of very small aspect ratio, for instance, 
this procedure would yield the lift distribution 

with a bending moment too low by 17 percent compared 
with the one of the exact lift distribution, 


cc,_A 
c — " 3 




whereas equation (26) yields 


cci_A 
c 4 


VKp)‘[ 1+4 (^)’] 


with a bending moment 8 percent too high. For very large 
aspect ratios, y D and y R both tend to the chord distribution 
c/c, so that the lift distributions given by equation (26) tend 
to the correct limiting value, 



(28) 


Therefore the accuracy of the results furnished by equation 
(26) may be expected to increase as the aspect ratio increases, 
whereas the accuracy of the other approximation is inde- 
pendent of aspect ratio. 

In this connection, it may be mentioned that “strip theory” 
consists in using equation (28) for all aspect ratios and, hence, 
is not very satisfactory for wings with medium and low aspect 
ratios. For instance, for the case discussed in the preceding 
paragraph, it furnishes a bending moment which is too low 
by 25 percent for a delta wing, and too high by 25 to 100 
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percent for tapered wings. (Tlie two figures pertain to taper 
ratios of ){ and 1, respectively.) 

Substituting the value for a from equation (27) into equa- 
tion (26) gives the following expression for cc,/c: 

{ [y B (y) -Ky s (y)+bKS(y- v )] V R (v) <h 

C 0 J-b/i 

The expression in braces in the integrand is the desired gen- 
eralized aerodynamic influence function and will be desig- 
nated by Q(y,v), so that 

<?(2M)=[TD(y)— K y R (y)+bK 8(y—y)] Tr(t?) (29) 

and 

rr, C L C bn 

G(y,v)«(y)d v (30) 

C o J-b/2 

In the limiting cases of wings of very low or very high 
aspect ratio, y D and y R approach a common value, say y. 
Thus, for wings of very low aspect ratio, 

C(y,v)=^ [T(y )+6 S(y—y)] 7 ( 17 ) (31a) 

where 5 is the Dirac delta (unit-impulse) function, and for 
wings of very high aspect ratio, 

G(y,y)=b 8(y~v) y(rj) (31b) 


which is the Green's function associated with strip theory. 
Thus, on wings of very high aspect ratio all the lift produced 
by the local angle of attack at a given station is carried in the 
immediate vicinity of that station, whereas on wings of very 
low aspect ratio much of the lift is carried elsewhere. This 
tendency for a given station to affect a great er portion of the 
wing as the span decreases is, of course, to be expected. 

The preceding analysis can be applied to the oscillatory 
case at a given reduced frequency Tc, as well as to the steady 
case. If the assumption of invariance of normalized distri- 
butions with time or frequency is made, as before, then 
y D (y), yn{y), and K are independent of frequency, so that 
equation (30) can be written as 

l(y;k) = ° La( ® qS j^GiV'VMridv (32a) 

where <9(y,i}) , defined as before, is independent of k. By 
applying a Fourier transformation to both sides of this equa- 
tion the following relation is then obtained for flight through 
continuously varying turbulence (cf. eq. (4)), as modified by 
the assumption stated in equation (5) : 

At(ti)dti|J* ^Q(y,ri)w(U(t— U),i))dri (32b) 

where l{y,t) is the lift per unit span at station y and time t, 
and where the function h L (i) is the one used previously. 


Tfie correlation function for this lift can then be written as 


'/'im ( r ) J" J h L (ti)hL(to)dtidti -p J ^ J ^ G(y,ili)G(y J T}2)ft tl , UU t (T+ti—t i ) t + drhdrjt 

=^ij^ J At(ti)Ai(<a)dtidti'/' tt . < (U (t— |— ^ j — tj), y) (33a) 

where the function 

'kv>,(P T ,y)=^jl J*_ ro Q(U,Vi ) G(y,T ]2 )f K .(-s l IU t T t + fa— ip ) s ) dijidria 


represents an effective correlation function, which when trans- 
formed into the equivalent power spectrum p Ut (u,y) can be 
used to obtain the power spectrum for l(y,t) and hence its 
mean-square value. Thus 

<Puv)(. 0} ) = fy U~) \4>(k)\ 2 <Pw,( C0 >y^ (33b) 

Before the calculation of <p 0f (w,y ) is discussed, the fimction 
^(Ur,y) will be defined in a somewhat more general form 
than in the preceding paragraph in order to anticipate future 
needs, namely 

'Pv t (UT,y u y s ) 

si 1 f G[y u vi) G(y 2 ,v2) 'k^ l U 2 T 2 + vi^dvi 

o J-h'iJ-bll 

If the assumption is now made that y R (y) and y D {y) are the 


same, as is the case for all unswept wings and for wings of 
very low or very high aspect ratio, they can both be identi- 
fied with the function y(y) used previously, so that 

G(y, v ) = [(l-K) y(y) + bK 8(y-y)} y(y) (34) 

and 

kw. (U r,yi ,y a ) = { ( 1 —K) Vir, (Ut)+K( 1—K) <JJr,yi) + 

K,- ( Ur,yi )] +i£Y.r(Vzy 2 T*+ (7/2— 2 /i) 2 ) } 7 (l/s) 

where 'P Ks {Ur) is the coiTelation function calculated pre- 
viously for the averaged vertical component of turbulence, 
and where 

\£'„ < .(27r,y)=|J' ^ 7 ( 77 ) ^„(V(7V+(y— >?)*) ch) 
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Similarly, then, 

‘Pa ( ( u >yi>yi ) = ( (i — 

<P*/(u,y i )]+K?p a (Q>,y a —y l ) }y(j/i)y(y 2 ) (35) 
where <p Ue ‘(u,y) is the Fourier transform of ^‘(Ur,?/), so that 

i r n 

‘Pu, e -(u,y)=frJ_^v(y)v>u,(u,\y—v \)dy (36) 


and thus represents an averaged form of the spectrum 
introduced previously. (See eqs. (16) and (19).) For uni- 
form spanwise loading and the point correlation function 
used previously, this function is 

11. 

<P<r e ‘ (") y)==^7" p { 2,k ,1 [Ki Q (ai) +Ki Q (a-i) — 

aiK 0 (ai) — a*Ko(aa)H- [ai*iT, (<m) +a^ s £Ti (aa)] } 

where 


and is shown in figure 4 for several values of jr^* 

The more restricted form of <p Ut required in equation (33b) , 
can now be obtained from the more general form given in 



oif/U 

Figure 4. — The power spectrum ?«,,•(<*, y). 


equation (35) by setting y 1 =y 1 =y in the latter, to yield 

<Pa t («) = [(1 — 2£)V» 4 («) + 2K{\ —K)<p Wt *(u,y ) + 

E?<p a (w)]y i (y) (37) 

MEAN-SQUARE BENDING AND PITCHING MOMENTS 

When the variation of the gust intensity along the span is 
taken into account, the mean-square lift L 1 is not an ade- 
quate index of the stresses in the wing, nor can the stresses 
be obtained from the mean-square lift distribution P{y)- 
Instead, the mean-square bending and twisting moments, as 
well as the mean-square vertical shear, must be calculated 
directly. 

In a manner analogous to that employed for the lift and 
rolling moment, each of these quantities can be expressed in 
terms of a certain influence function which, by means of 
the reciprocity theorem, can be related to a lift distribution 
on the wing in reverse flow. For instance, for the root 
bending moment the desired lift distribution is the one for 
an angle of attack which is zero on one wing and proportional 
to y on the other, as may be seen from the fact that the bend- 
ing moment can be expressed as 

n 

[l(y)]ai,)j(3l)dy 

-M 5 

where 

f(y)= 0 (y< 0) 

f(v) =y (y^o) 

so that, according to the reciprocity theory for unswept 
wings in steady or indicial flow, 

[l{y)]f(.vMy) ( h 

-W1 


Similarly, for the root shear the required lift distribution 
corresponds to an angle of attack which is zero on one -wing 
and uniform on the other. 

Such lift distributions can be calculated readily. For 
instance, for subsonic flow and unswept wings they can be 
obtained from those given in reference 10; the lift distribu- 
tion for the root bending moment is one-half the sum of a 
linear symmetric and linear antisymmetric lift distribution, 
and the lift distribution for the root shear is one-half tho 
sum of a lift distribution due to a uniform angle of attack 
and a lift distribution due to deflection of a full-span aileron. 

However, in some cases such calculations may be time- 
consuming, and an alternative approach may be desirable. 
One such approach consists in synthesizing the desired influ- 
ence function from the generalized lift-influence (Green’s) 


function discussed in the preceding section. In order to 
illustrate the use of this approach, it is adopted in this section. 

The bending moment at any station y(o^y^-^ and at 
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any time t can be obtained from the lift distribution considered in the preceding section as 


r V» 

M B (y,t)= J (y'—y)l(y',t)dy' 

i rw c » i r>>n 

=j J <y— 2/)dj/J_ v OW,v)iDm-tMdn 

= l j- ^&{y,-q)w(JJ {t—t x ),r,)d.T) (38) 


where the influence fmiction for the bending moment is 

=jJ (y'—y)G(y',v)dy' (3i)a) 

so that, upon introducing the previously used function for G(y,v) 

M(y, v )=[(l-IOM 1 (y)+KM 1 (y, v )}v( v ) (39b) 

where 

i r 6/s 

Mi (y) J (y'—y)y(y') dy 1 

Mt(y,y)=y—y (v>y) 

M 1 (y,T ])= 0 (ij^-y) 

and, hence, for the root bending moment, 

Mi(0)=l£ n y'y(y')dy' 

Mt(0,rj)=i} (y>0) 

M,( 0,„)=0 (ujgO) 


Hence, the correlation function for the root bending moment is 

X a> pm 1 /'ft/a /’&/! 

h L (ti)h L (.Widh± {a-K)*M 1 *(y)+K(l-K)M 1 mM 2 (0,vi) + 

-CD J — CD 0 J -k/2 

M a (0,y^)]+E}M 3 (0 ) y^M a (0 1 y i ) }y{ni)y{v^'i'wblU 1 {r+ti—U) iJ r{vi~yi) v )dyidyi (40) 


The approach used in the preceding sections can now be used 
to obtain the power spectrum of the root bending moment 
by evaluating the inner pair of integrals of equation (40), 
taking the Fourier transform of the result with respect to 
r, and multiplying the power spectrum obtained in this 


manner by 


(Sg0W)l* or by 


approaches indicated for the lift. 


using one of the alternative 


The spectrum corresponding to the inner pair of integrals 
is, for K~ 1 and uniform loading, 


(32/3'+0' 3 )#i (j)]+fc' s [f*Ki 0 (|')+(64-6 J 3 ,i )- 

W*K 0 (£)- (32/3 f-pT^Kt (I)] y 
where /3 / =^/i+p3 ; as before. 


The integral of this spectrum is 

^.”( 0) =^[(- 24 +f)+e-^ ( 24+ 12 d+ 3 d 2) ] 

Examination of ^„ t /,( 0) indicates that, although the moan- 
square bending moment tends to decrease as the span in- 
creases, it decreases less rapidly than the mean-squaro lift, 
with the result that the effective lateral center of pressure 
moves outboard. Quantitatively, the distance from the 
plane of symmetry to the effective lateral center of pressure 
can be defined as the square root of the ratio of’ the mean- 
square root bending moment to the mean square .of the lift 
on one wing, that is, of the root shear. Although those mean 
squares have not been calculated, the square root of the 
ratio of ^„ e #/(0) to the corresponding value for the root shear 
increases by 15.5 percent as /S increases from 0 to infinity, 
with much of the increase realized at fairly small values of fi. 

For a swept wing the variation of the gust intensity along 
the span results in a pitching moment which must be taken 
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into account in calculations of the dynamic response of the 
airplane to continuous turbulence. This pitching moment 
can be obtained in substantially the same manner as the 
bending moment. Thus, if y is the station of the mean aero- 
dynamic chord, 

J in 

(y—\v\)i(y,t)dy 

-i/i 

i p® ps/i 

=tanA^J A £ (ti)dtiJ ^M s (y)w(U(t—ti),y)dT] (41) 
where 

i rm 

M 3 (v ) J (y— |y| ) Q(y,y)dy 

so that, with the previously used approximation to the 
Green’s function, 

' M 3 („) = { (l -K) \y-2MM\ +K®-\n |) Mu) 

The correlation function, spectrum, and mean-square value 
of the pitching moment can then be obtained in the manner 
used in the preceding sections. 

WING-TAIL CORRELATION EFFECTS 

The tail strikes a given gust some time after the wing does; 
as a result, a pitching moment arises which does not exist in 
steady or quasi-steady motion, nor if the airplane' is very 
small, because then the time lag is insignificant. This pitch- 
ing moment can, for the purpose of the present paper, be 
analyzed either in terms of the correlation between the gusts 
at the wing and those at the tail or, if a time-lag term is 
included in the indicial-response function of the tail, in terms 
of the correlation between the wing and tail response func- 
tions. The first point of view serves to exhibit the effect 
under consideration more clearly and is adopted first; the 
second is more convenient and is adopted in the subsequent 
parts of this report. 

In the somewhat artificial case of a small wing and tail 
separated by a relatively large distance, only the distribution 
of turbulence along a line (the flight path) rather than in a 
portion of a plane is needed. The pitching moment due to 
the vertical component of atmospheric turbulence can then 
be written as 

M(t)= f h a (ti)w(U(t—t 1 ))dt 1 + f htitijwixt+Uit—ti^dti 

J —CO J —to 

where x , is the tail length, and where h„(t) and h,(t) are the 
pitching-moment responses to indicia! gusts hitting the 
wing and tail, respectively, at <=0; both may include un- 
steady-lift effects, and, if downwash effects are to be con- 
sidered, h w (t) should include the contribution to the pitching 
moment of the tail lift caused by the downwash at the tail 
associated with the lift on the wing which results from the 
iudicial gust. 

The spectrum of this moment can then be written as 
<pjr(«) =jjB»(&>) +H,(' m)|*— 2R-j^l— e^ 1 '^S' a ,*(co)B' I (a))^-Jp Kr (&j) 


where the symbols R{ } and * designate, respectively, the real 
part and the complex conjugate of a complex number. 
Obviously, when x, approaches 0 the second term in the 
bracket in this equation vanishes, so that the first term repre- 
sents the perfect-correlation effect, and the second represents 
the correction for imperfect correlation. 

In order to furnish an estimate of the magnitude of the 
effects under consideration, some calculations have been 
made on the basis of the assumption that a real coefficient n 
exists such that 


H w (u)=nH t (cS) 

which implies that the attenuation with frequency of the 
contributions of the wing lift and tail lift, respectively, to 
the pitching moment is the same. The ratio ji is —1 for 
neutral stability, and ;C> — 1 for stable flight; it is positive 
when the aerodynamic-center location (tail off) is behind 
the center of gravity, so that positive values of /i are not 
likely to be incurred with normal configurations and flight, 
conditions. 

For this case, 

’hr (<*>) = | H,(w) I 1 jjl-f m) s — %n(l— cos <Pw (<■>) 


and 


M} 


(M 


(tail alone) 


.(l+„)>-2 ,*(fi.f|) 


where 



|fit(«)|Vu>(«) 1 

'o 

(i-cos^ 

| do> 

l~ 

!#<(") \ V»(<*>)da> 

0 


The function has been calculated for several values of its- 
arguments by means of the lift-attenuation function given 
in equation (7) and the point spectrum used in the pre- 
ceding sections, and is shown in figure 5. Also shown is the 
ratio of the mean-square moments as a function of the 
factor n which in figure 5 is referred to as MJM t . The effect, 
of imperfect correlation is seen to be very large as the- con- 
dition of neutral stability is approached; the entire pitching 
moment is then the result of instantaneous differences in 
gust intensities at the wing and tail. 

In general, however, the tail length and the span are of the 
same order of magnitude, so that an analysis of the effect of 
imperfect correlation between the wing and tail must take 
into account the averaging effect of the wing span. The 
pitching moment at any instant is then 




7(y) w(U(t—ti),y)dy+ 



h t (ti)w(x t +U(t—U)fi)dti 


(42a) 
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Figure 5. — Effect of tail length x t on the mean-square pitching 
moment. 

Hence, the power spectrum of this moment is 
Ojf (w) = |H„(o>) I (") + 1 j St (<o) I Vw(“) + 

2R-j^e v (42b) 

where ^.(w) is the averaged spectrum of equations (11) 
and (13), and <p Wt -(w,0) is the spectrum of equation (36) for 
y=0. (It should be noted that h,(i) is the response to an 
indicial response which strikes the tail at <=0; if it were the 
indicial response to a gust which strikes the wing at f=0, 

the factor e ^ in the preceding equation would not be 
required.) 

DYNAMICS OF THE RIGID AIRPLANE 
In this part of the report the motions of a rigid airplane 
subjected to atmospheric disturbances are considered. This 
part thus applies directly to those airplanes which fly at 
relatively low speeds and do not experience any significant 
structural deformations, and also serves as a preliminary to 
the treatment of the flexible airplane in the later parte. 

The assumption is made that the motions are small enough 
to permit the use of linear approximations to the resulting 


aerodynamic forces and to permit the linear superposition 
of these forces. The longitudinal degrees of freedom 
(pitching, vertical, and horizontal motion) and the lateral 
degrees of freedom (yaw, sideslip, and roll) can thoreforo 
be considered separately. The first two sections of this 
part are concerned with the longitudinal motions of a small 
airplane. The material presented here is thus a generaliza- 
tion of the single-degree-of-freedom analysis presented in 
reference 2. Mathematically, the dynamic system is now 
described by a set of simultaneous ordinary differential 
equations rather than a single one, but the problem of calcu- 
lating the pertinent transfer functions is still one of simple 
algebra. The extension of these results to the largo airplane 
is effected in the third section, using the techniques developed 
in the first part of this report. The lateral motions of a 
large airplane are considered briefly in the fourth section, 
and the combination of the stresses due to longitudinal 
and lateral motion of a large airplane is discussed in the 
fifth section. 

EQUATIONS OF LONGITUDINAL MOTION 

The equations of motion of an airplane can be expressed 
in several coordinate systems. The system which is gen- 
erally the most convenient one for analyses of airplane 
stability and is generally referred to as stability axes, con- 
sists of body-centered axes which are normal and parallel 
to the relative air velocity and rotate with the airplano as it 
pitches or yaws. (See ref. 11, for instance.) The aero- 
dynamic forces related to this axis system can bo measurod 
more readily in wind tunnels than those related to other 
axis systems. In view of the very close relation of a stability 
analysis to the problem considered here, these axes will be 
used in this part of the paper, but in the analysis of the flexible 
airplane in the subsequent parts of the report space-centered 
axes will be used, because they are slightly more convenient 
for that purpose. 

The airplane will bo considered to be in steady level flight 
prior to disturbance. The motions studied will be the de- 
viations from their mean values; for instance, the angle 0 con- 
sidered here will be the difference between the disturbed and 
the initial value of the angle of pitch. Hence, the motions 
and forces calculated by the method indicated here must bo 
added to their mean values to obtain the total motions and 
forces. 

Inasmuch as, for the purpose of a statistical analysis, tho 
dynamic characteristics of the airplane are represented most 
conveniently by its transfer functions, attention will bo 
confined in this section to sinusoidal gusts and motions. 

For this case the linearized equations of longitudinal 
motion can be written as follows (see eqs. 11-193 of ref. 1 1 , 
for instance) : 


Ka — Z a 

-Z u 

— iuU 


"to," 

— X a 

ico — X u 

9 

-> 

Up - 

_ — VaM^— 

-M u —M u 

— CO 2 — icoM, _ 


. e . 


=4>(k) 


YZ« 

X u X u 
[_iW„ M u J 


O 


(43) 
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The stability derivatives which appear in these equations 
are defined in table 1 in terms of conventional aerodynamic 
coefficients, and the numerical values are given for the ex- 
ample used in reference 11. (The value of the mean chord 
is not given in ref. 11 but it is assumed herein to be 10 feet 
on the basis of other information given in ref. 11.) 

In analyses of the stability of a rigid airplane the quasi- 
steady approximation to unsteady-lift effects is usually made, 
in which, in effect, the forces corresponding to a steady 
attitude, to constant disturbance velocities, and to constant 
accelerations are considered. This approximation is justified 
because the motions of concern are generally sufficiently 
slow. For the same reason this approximation can also be 
made in analyzing the response of an airplane to atmospheric 
turbulence. 

However, in this problem another type of unsteady-lift 
effect occurs, namely, that related to the forces directly 
attributable to the turbulence. This effect is here taken 
into account by multiplying the quasi-steady values of the 
forces due to gusts on the right side of equation (43) by the 
attenuation function This procedure implies the 

assumption that the airplane is small relative to the scale 
of turbulence, inasmuch as no averaging effects have been 
taken into account; these effects will be discussed presently. 
Also, this attenuation function is strictly applicable only to 
the normal forces. 

The unsteady effects on the drag are not known because 
of the relatively complicated nature of the mechanism which 
gives rise to drag. If, however, the assumption is made that 
upon entry into a sharp-edge gust the drag rises linearly and 
attains its steady-state value in the time required to travel 


TABLE 1.— DEFINITIONS AND NUMERICAL VALUES OF 
STABILITY DERIVATIVES 


Derivative 

Definition 

Numerical value for 
the example 


~&( Cl ° +Cd ) 

— 1. 430 

Zu 

-Sf?( c ‘.+ c -> 

-0. 0955 

A r „ 


0. 0016 

Xu 

~tS( C d - +Cd ^ 

-0. 0097 

M. 

qSc r 

-0. 0235 

Mi 

qs& r 

2/y„U J " 

-0. 0013 

M u 

B Tt C - +C -> 

0 

M q 

qSc* p 
2 I V ,U "■> 

— 1. 920 


W, lb 30, 500 

U, fps — 660 

c, ft 10 

Altitude, ft 20, 000 


526597 — 00 12 
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N chord lengths, the drag equivalent of |<£(&)| 2 is the function 

1— cos 22V7: 

2W 

which, for N equal to about 5 or 6, agrees fairly well with 
|<£(&)| 2 in the region of main interest (&>1). 

The unsteady moment is also difficult to predict because of 
the paucity of knowledge concerning unsteady downwash 
effects for wings of finite span. However, inasmuch as the 
wing lift contributes part of the moment and, through the 
mechanism of downwash, determines to a large extent the 
moment contributed by the tail, the use of the lift attenua- 
tion function for the moment appears reasonable for a first 
approximation, and the use of the same function for the lift, 
drag, and moment facilitates the analysis. 

For a more refined analysis, the time lag between the 
instants at which the gust hits the wing and the tail must be 
taken into account, not only for large airplanes, but even for 
small airplanes if phugoid motions are important. (See ref. 
12.) Also, the lag in downwash should be taken into 
account for large airplanes and possibly also for small air- 
planes in some cases. One way of achieving this result is 
indicated in the next part of the present paper; another is 
discussed in reference 12. 

In equation (43) the unknown quantities w p and u v are 
the normal and axial components of the disturbance velocities 
of the airplane relative to the free stream. Inasmuch as the 
coordinate axes rotate during the motion, the time deriva- 
tives of these quantities do not represent the actual airplane 
accelerations, which are required in analyses of the loads 
experienced by the airplane and the degree of passenger 
discomfort. H the deviations from a mean flight path are 
assumed to be small, the accelerations normal and parallel 
to the chord or longitudinal axis of the airplane are sub- 
stantially the same as the absolute vertical and horizontal 
accelerations z and x, which can be obtained from the 
relations 

z=-w P -\-U6 

x=—u P —gd 

The transfer functions for these quantities can then be ob- 
tained by introducing these relations into equation (43). 

In studies of the longitudinal stability of airplanes, equa- 
tion (43) is rarely solved in the form given here. It is 
usually reduced to two equations with two unknowns, 
either Up and 6 (the phugoid case) or w P and 6 (the short- 
period case), the short-period case being usually the one of 
primary interest. The part of the turbulent energy con- 
tained in the frequency range near the phugoid frequency is 
relatively small, so that the phugoid case has no significance 
for the analysis of loads and accelerations resulting from 
atmospheric turbulence. Hence, the short-period case, 
which ignores the phugoid oscillations, furnishes an ex- 
cellent approximation to the loads and accelerations asso- 
ciated with the longitudinal motions of an airplane in turbu- 
lent air. However, another two-degree-of-freedom case, 
the one involving w T and u P , is useful in certain studies of 
the effects related to the interaction of horizontal and vertical 
components of turbulence. 
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Both of these two-degree eases can be reduced to the single- 
degree-of-freedom case involving only z (or w P ). For air- 
planes which have a large moment of inertia in pitch this 
simple case furnishes a good approximation. It has been 
studied in reference 2, where substantially the same approxi- 
mations to the unsteady-lift effects were made as are made 
here, except that in reference 2 apparent-mass effects were 
included. (These effects are not included in the stability 
derivatives used in equation (43) because they are usually 
small — less than 1 percent of the mass of the airplane — and 
are different for each degree of freedom. However, if 
desired, the apparent mass pertaining to a given degree of 
freedom can easily be added to the airplane mass in calcu- 
lating the stability derivatives.) However, this approxima- 
tion is more nearly valid for calculating peak loads pursuant 
to an entry into a sharp-edge gust than for calculating the 
response due to random turbulence. Consequently, in the 
following sections, attention will be confined to the short- 
period case, although the analysis is equally applicable to 
the other case and easily extended to the case of three degrees 
of freedom. 


SOLUTION OF THE EQUATIONS OF LONGITUDINAL MOTION 


Transfer functions. — If the degree of freedom pertaining 
to x (or ,u p ) is ignored, the solution of equation (43) can be 
written as 


(w) Hi (to) 

X'eJ _H?(«)H?(«) 



(44) 


where the transfer function H^(w) is defined by 




— Co 


(45) 


where, in turn, the coefficients are defined in terms of the 
stability derivatives (see table 1) by 


except that Z„, ZJ , M m and MJ are replaced by Z„, Z u ', 
M u , and M u r . 

With these transfer functions the mean-square values of 
2 and 6 can be calculated from the spectra of w and u by 
using equation (1), provided that the simultaneous action 
of w and u is taken into account. In order to analyze this 
effect the vertical acceleration z will be considered, but the 
analysis will be applicable to 8 or any other characteristic of 
the airplane which responds to w and u. Furthermore, the 
transfer functions need not be those considered in equation 
(44), but can be those calculated for the three-degree-of- 
freedom system or for a flexible airplane. 

Normal-acceleration spectrum and mean-square value. — 
For the present purpose the indicial-response functions 
hi it) and hi it), which are the Fourier transforms of the 
transfer functions (see eq. (3a)), are more convenient. In 
terms of these indicial-response functions, 

2(f)= f h%{ti) , w(U (t—u)) du+ [ tii(ti) u(JJ(t— u)) (Hi 

• J — ' CO 1/ CD 

Then, if wit) and u{t) are stationary in a statistical sense, the 
correlation function for z(t) can be written as 

<Pi(r)= J J [A'i(ii) h'liti) ^'»(D r (T-|-0~0))+ 


hiiti) h;it 2 ) f u (U{T+U—t^) + 
h'lit i) h'liU) &™(Z7 (r+0— 0)) + 

hiiU) hiiti) twuiUi— T+ti— 0 ))] dti dtt 

(46) 

where ^w(Ht) is a cross correlation of w and u defined by 
'pwuiUr)=wix J rUt) uix+Uit+r)) 


A 0 =l 

Bo= — (Zuj-f-iWj-f- ZJM &) 

Ca=hl^Z w — Uhl u , 

The transfer function 2?j(co) 
terms of the coefficients 


Al=-Z w r 
Bl=Z a ' (UMu+M,) 

01= —U(MJZ a —MJZJ) 
can., be defined similarly in 


Now, if the turbulence is isotropic, the mutually perpen- 
dicular velocity components u and w at points in the .XT- 
plane are statistically independent, and their cross correlation 
is zero. Therefore, the two terms in equation (46) involving 
'i'wuiUr) vanish, and the power spectrum of z is 

« («)“ | Bf («) | W»)+ 1 Hi («) IV» (47) 


Al=M a '+ZJM- w 

B%=-(hl^Z w -MJZ^) 

a = o 

In these equations a distinction has been made between the 
values of Z a and M w which occur on the right side of equation 
(43) and are here designated by a prime mark, and those on 
the left side of that equation. The primed derivatives 
pertain to the lift and moment directly due to gusts, whereas 
the unprimed derivatives pertain to the lift and moment due 
to airplane motion; the reason for this distinction is discussed 
in a later section. Furthermore, the coefficients A%, Bl, C%, 
Ag, Bg, and C| are the same as the coefficients A ”, J5|, . . . 


so that, generally speaking, the power spectrum of a response 
which depends on both the horizontal and the vortical com- 
ponent of turbulence is simply the sum of the power spectra 
of the two contributions, provided the turbulence is isotropic. 
(This statement can be shown to be true even if the distribu- 
tion of the gusts over the span is taken into account.) 

For the short-period two-degree-of-freedom case, then, the 
contribution due to w is 


2 — 



^-[2Aici-(Biywi-(oiy 

A-pCo-BaW+Co 1 


<Pa( U)d 


03 


and the contribution due to u can be obtained from the same 
expression, but with the subscript and superscript w replaced 
by u. However, the ratio of the two contributions is in the 
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order of 4a 3 , where a is the trim angle of attack in radians, 
measured from the zero-lift conditions. Consequently, ex- 
cept at very high lift coefficients, such as those used in land- 
ing, the contribution due to u is usually negligible compared 
with the one due to w and is disregarded in subsequent 
sections. 

By means of the lift-attenuation function given in equa- 
tion (7) and the point spectrum used in the first part of the 
report, the integral for jj 3 has been evaluated (by using the 
technique of partial fractions for the integrand) for the ex- 
ample of reference 11. (The lift-attenuation function of 
equation (7) has been used despite its shortcomings in order 
to facilitate the analytic integration of the spectra.) The 
results are shown in figure 6, as are the results calculated 
similarly for the three-degree-of-freedom case, the other 
two-degree-of-freedom case (horizontal and vertical motion; 
referred to in the figure as the zero-pitch approximation), 
and the single-degree-of-freedom case (vertical motion) . An examination of this figure indicates that, at least for this air- 

plane, inclusion of horizontal motions does not affect the mean-square normal-acceleration response to any significant 
extent. These calculations pertain to a “small” airplane, and no distinction has been made between the primed and unprimed 
values of Z u and M a . 

The preceding treatment of the short-period case has the advantage of using readily available information concerning 
the characteristics of any given airplane. For the purpose of trend studies a dimensionless form of the transfer func- 
tions is preferable. The preceding equation for s 3 can be written in dimensionless form (the contribution of horizontal 
gusts again being neglected) as 



L *, feet 

Figure 6. — Effect of scale of turbulence on the mean-square 
acceleration of the example airplane. 


and, similarly, 


ZlY^i fVrMi* ^ +4 ( y U<p w {cS) ■ 
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(2 gfcY~ \gcf2j ZyLy k) +h J J 0 '^ ( )l k*-2 (fc 0 3 -^)P+(fc 0 3 + y 3 ) 3 L*v? 


(48a) 


(48b). 


where k is the mass parameter 


8m 


K =C LaP Sc 

v is the dimensionless damping coefficient 

_ c log, 2 

v 2U Ta 

Ty u is the time to damp to one-half amplitude, which is 
given by 

l °M^=-±(Z w +M t +UMJ 

so that 
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and k 0 is the dimensionless frequency of the short-period 
oscillations, 
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Thus, for this two-degree-of-freedom case the dimension- 
less mean-square responses are functions of only two addi- 
tional parameters, which are dimensionless forms of the main 
characteristics of the short-period case (the short-period fre- 
quency and the time to damp to one-half amplitude), beyond 
those encountered in the single-degree-of-freedom case, 
namely, the mass parameter k and the scale parameter 
_L* 

S ~c/2 

Bending-moment spectrum and mean-square value. — For 
a small airplane the instantaneous bending moments at vari- 
ous points on the span are proportional to the instantaneous 
normal acceleration. For instance, the root bending moment 
can be written as 


A ,0 — Tfl w — _ 

B= ~2 V ~~2 yz 


(49a) 
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where L a and m„ are, respectively, the lift on and the mass 
of the wing, and where y and y are the lateral distances to 
the center of pressure of the lift on one wing and the center 
of gravity of the mass of one wing. However, L a is propor- 
tional to the lift L on the entire airplane, which in turn is 
proportional to the normal acceleration, so that 



and 

"*-[! (%Hsr# < 48b > 

Hence, the spectrum and mean-square value of M B are pro- 
portional to the spectrum and mean-square value of 2 , re- 
spectively, the constant of proportionality being the square 
of the quantity in the brackets of equation (49b) . 

SPECIAL PBOBLEMS BELATED TO THE LONGITUDINAL MOTION OF LARGE 

AIRPLANES 

Single-degree-of-freedom case. — In the preceding sec- 
tions the airplane has been assumed to be small in the sense 
of this report, and neither the instantaneous lateral variation 
of the intensity of turbulence nor the difference between 
instantaneous intensities at the wing and the tail has been 
taken into account. In this section this restriction is 
removed by introducing the aerodynamic forces calculated 
in the first part of this report into the dynamic analysis of 
the preceding sections. The arguments advanced in the 
preceding section for ignoring horizontal gusts and horizontal 
motions are equally valid for the large airplane; therefore 
these gusts and motions will not be considered here. 

For. the single-degree-of-freedom case involving only 
vertical motion the required modification for the normal 
acceleration is very simple. For this case the transfer 
function is 

h?( (so) 

wJ ^ to x 


where ZJ , attenuated by is the stability derivative 

for vertical gusts and, hence, represents the lift per unit 
gust intensity. Therefore, if the result for the lift calculated 
in the first part of this report is used, the mean-square 
normal acceleration becomes 


(The function <p w (co) is defined in the first part of this 
report.) This expression differs from the result obtained 
in reference 2 only in that <p a (oj) is here replaced by ^(u). 

However, even for the single-degree case the calculation 
of the bending moment now becomes a considerably more 
complicated problem, because the lateral centers of pressure 
of the lifts due to the motion of the airplane and directly 
due to turbulence no longer coincide. Equation (49a) now 
becomes 


M b = 


y L„ Oz a iS 
2 L U 


z 

U 



1 /*■» rw 

6 J-„ h L ^ dy (6i) 


where 0 La is the lift-curve slope for the entire airplane, 
h L (t) is the response function used in equation (38), and 
y&i(y) is the function i&(i/,»?) used in equation (38), with y— 0. 
In the first two terms on the right side of equation (61), 
2 and 2 can be related to w by means of h“ (the Fourier 
transform of the function presented in eq. (50)) and the 
lift directly due to turbulence obtained in the first part of 
this report, so that these two terms can be written in the 
form 

1 f “ C6/2 

6 J dti J w{ ^ d y 

and, hence, equation (51) becomes 

J dri [Al(£i) Ta r(j/) 

hM(ti)y{y)]w(U{t— ti),y)dy (62) 

The required influence function for the bending moment 
is, thus, 

hx B (y, t ) [h L (t) V M (y) +h M (t) 7(y)] (53) 

Although the two terms of this function are products of a 
time-dependent and a space-dependent constituent, as in 
equation (5), their sum cannot be split up in this manner. 
Hence, the techniques used in the first part of this paper 
are not directly applicable. The general treatment of prob- 
lems involving influence functions for which the assumption 
of equation (5) is mot valid will be considered in the last 
part of this paper. However, in this section a special 
technique will be used that applies to cases for which the 
influence function can be expressed as a sum of several terms 
(two in this case), each of which can be expressed as a 
product of two functions, which depend, respectively, on 
time alone and distance along the span alone. (See also 
ref. 13.) This approach is more convenient than the general 
approach of the last part of this report when the number of 
terms is two or, possibly, three. (Although in this section 
only two terms will be considered, the generalization to 
three or more terms is straightforward.) When the number 
of terms is greater than three, the general approach becomes 
more covenient. 

. If the value of M s given by equation (52) for t+r is 
multiplied by the value for t and the result is averaged, and 
if the Fourier transform of the resultant correlation function 
is then taken, the following expression is obtained for the 
spectrum of the bending moment: 

<Pm b (o>) = ISx(co) + \S M (u) I («) + 

2-R{H L (oS)H M *( u )} ip „' i (cS) (64) 

where the symbols R{ } and * designate, respectively, 
the real part and the complex conjugate of a complex num- 
ber. In this equation H M {u) and H L (w) are the Fourier 
transforms of h M (t) and h L (t), v c< (w) is the previously defined 
averaged spectrum for the lift, and <p v ^ and <p v ^ are averaged 
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spectra obtained in a similar manner but with the value of 
T given by equation (12) replaced by 

2 rm- n 

r, (’ 7 ) == 6 J-mz y ^ 7if ,(y+v)dy 

2 f<W- 1 " ^ 5 ) 

r 3 (>;)=T [yir,(y)vM,(y+v)+yif a (y)'Yv a (y+v)]dy 

u J -6/2 

J 

where the subscripts s and a refer to the symmetric and anti- 
symmetric parts of y M , respectively. 

The contribution of the antisymmetric part of y^ to r 2 
and, hence, to the spectrum for the bending moment stems 
basically from the asymmetry of the instantaneous distribu- 
tion of gust intensity over the span. This asymmetry gives 
rise to a rolling moment (which was considered in the first 
port of this report) and, hence, to rolling motions, which 
contribute additional bending moments due to the aero- 
dynamic and inertia loads associated with these motions. 
If, for the purpose of calculating the bending moment due 
to symmetric flight through turbulent air alone, these mo- 
tions are disregarded, then the contribution of y M to r 2 
should be disregarded as well. The problem of combined 
symmetric and antisymmetric motion will be considered 
presently. 

As pointed out in reference 13, the mean-square bending 
moment calculated in this manner may be smaller or greater 
than the value calculated by ignoring spanwise variations 
in gust intensity. If the mass of the airplane is almost 
entirely contained in the fuselage, the decrease in the lift 
which results from taking these variations into account 
causes a decrease in bending moment. However, if most 
of the mass is in the wing, the net bending moments (aero- 
dynamic less inertia) for a uniform spanwise gust are very 
small, and the effect of taking spanwise variations of gust 
intensity into account is to increase the mean-square bending 
moments. 

Two-degree-of-freedom case. — For the two-degree-of- 
freedom (short-period) case the analysis given for the normal 
acceleration in the preceding section can be extended as 
follows: As indicated in equation (44) the transfer function 
for z is now 



1 <o l (UMp, -f-Aig) i<o -\-UM w , 


1 . U(-Z„) 

I y (J; 


r I y MJ <f>{k) 

0 _ 


—ZJ) 4>(k) + 
(56) 


(The following analysis can be applied equally well to § by 
using Hf (w) instead of Hf(co).) 

If the expressions inside the brackets of equation (56) are 
designated, respectively, by Hi(<o) and and their 

Fourier transforms by h t (t) and h- 2 (t), then, as a result of the 
definitions of ZJ and MJ , 


z(t)= f hi{ti)L(t-ti)dti-\-[ hi(t{) dti (57) 

J — CO 4 / —(O 

where L(t) and M(t) are the instantaneous lift and pitching 
moment due to the vertical component of atmospheric 


turbulence, which have been obtained in the first part of 
this report. The calculation of pj( «) or of 2* thus requires 
not only the spectra of L and M, the calculation of which 
has been discussed, but also the cross spectrum of L and M, 
which has to be calculated directly from equations (4) 
and (42a) . The result is 

<P ~2 (“) = 1-5) (co) I 1 |H(co) |V»,(«) + |5i (w) |* [|S^„(co) | a <Pu, e (oi) + 

!5^&>)|V»(«)+2R-^e ^ 77, (co) H w *(<o) v’ lV (“ ; 0)J+ 

2 [R{f?i(w) £T 2 » H(«) F B *(u)}v> ff# (u) + 

R { Bi (co) H„*( co) E(u) H t *( <o ) } <p W( , (co, 0) ] (58) 

where the first two terms represent the contributions of the 
spectra of L and M, respectively (see eqs. (14) and (42b)), 
and the third represents the contribution of the cross spec- 
trum of L and M. For the present purpose the functions 
H((o), H k (o>), and H t (co) can be expressed as 

H(<o)=m( — Z w ) <f>(fc) 

H w (co)=m Ax a (—Z w ) <j>(k ) 

77, (co) = to Ax a (— ZJ] <t>(k ) 

where Ax a is the distance from the aerodynamic center (tail 
off) to the airplane center of gravity: This definition of the 
contributions of the wing and tail to the pitching moment is 
based on the considerations that the direct contribution 
of the wing can be estimated with good accuracy and the 
total pitching moment is likely to be known from experi- 
ments, so that the contribution of the tail (which includes 
the effect of the wing lift on the downwash at the tail) can 
be determined as the difference of the two. The functions 
<Pu, e (co) and ?>»,•(«, 0) have been defined in the first part of 
this report. In view of the fact that the function <f>(k) con- 
tained in some of the terms of equation (58) always appears 
in terms multiplied by others which contain 4>*(k), only the 
absolute square of this function is required, as before. 

SPECIAL PROBLEMS RELATED TO THE LATERAL MOTION OF LARGE 

AIRPLANES 

The equations of motion in the lateral degrees of freedom 
(roll, yaw, sideslip) have the same form and can be solved 
in the same way as the equations for the longitudinal motion. 
(See pp. HI-53 to HI-67 of ref. 11.) Again it is convenient 
to cast the problem in the form used in a stability analysis 
in order to take advantage of the results of such an analysis. 
For a small airplane it is necessary only to replace the terms 
due to rudder deflection by corresponding terms involving 
side gusts, namely, 


Yl B 8 R bjY,£<f,'(k) 
Ls r 8 r by Jj?jj4>'(k) 
N Sr 8 r by N^4>'(k) 
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where is a suitable side-force attenuation function for 

side gusts and the notation of reference 11 is used for the 
other terms. The terms corresponding to aileron deflection 
can be disregarded for the small airplane. In the lateral 
degrees of freedom the small airplane thus reacts only to 
side gusts. On the other hand, the large airplane also reacts 
in the lateral degrees of freedom to vertical gusts, through 
the rolling moment calculated in the first part of the present 
report. If this rolling moment is to be included, it replaces 
the term L Sa <5 a used in reference 11. 

Instead of treating all three degrees of freedom simul- 
taneously, in stability analyses two one-degree-of-freedom 
cases are often considered, namely, the one of sideslip alone, 
with angle of yaw equal and opposite to angle of sideslip 
(the Dutch roll case), and the one of rolling alone. 

The Dutch roll case may be used for gust-load purposes 
in connection with yawing and sideslipping motion due to 
lateral gusts, provided the phase of the motion is not im- 
portant. (As may be noted from the preceding sections, 
the phase of a transfer function is important only in terms 
involving cross spectra.) Also, the Dutch roll case may 
prove useful in calculating the vertical-tail loads resulting 
from flight through turbulent air, particularly if the human 
or automatic pilot holds the wings substantially level by 
means of the ailerons. However, in general it does not ap- 
pear to be as satisfactory an approximation as the one- 
and two-degree-of-freedom approximations for longitudinal 
motion. 

For rolling motions due to rolling moments, the single- 
degree case of rolling alone appears to furnish a very good 
approximation. Although the rolling motion causes yawing 
and sideslipping motions, these motions do not appear to 
reflect on the rolling motion. Thus, the rolling motion 
which results from the rolling moment can probably be 
calculated fairly accurately without regard to the other 
lateral degrees of freedom. Furthermore, within the assump- 
tion of small motions, the stresses associated with these 
other lateral degrees of freedom do not generally contribute 
appreciably to those associated with the longitudinal degrees 
of freedom in the parts of the structure for which the latter 
are critical, such as the wing (although they may be critical 
for other parts of the structure, such as the vertical tail). 
Therefore these degrees of freedom (yaw and sideslip) will 
be ignored in the treatment of the large flexible airplane in 
the last part of this paper. However, if chordwise bending 
effects (deformations parallel to the chord) are important, 
as they may be in some cases at speeds close to the flutter 
speed, these other lateral degrees of freedom may have to be 
included in the analysis. 

For the large airplane, which responds in the lateral 
degrees of freedom both to vertical gusts and to side gusts, 
the superposition of the resulting responses, such as stresses, 
may be effected in the way indicated for the' interaction of 
horizontal and vertical gusts. If the turbulence is isotropic, 
the vertical and lateral gusts are statistically independent 
for points in the AY-plane, so that the spectrum of a given 
response is equal to the sum of the spectrum of that part of 
the given response which is due to vertical gusts and the 


spectrum of that part of the response which is due to lateral 
gusts. 

COMBINATION OF THE RESULTS OBTAINED FROM ANALYSES OF THE 
LONGITUDINAL AND LATERAL DEGREES OF FREEDOM 

The instantaneous wing stresses depend both on the 
motions in the longitudinal degrees of freedom (primarily 
vertical motion and pitching) and on those in the lateral 
degrees of freedom (primarily rolling) . The purpose of this 
section is to indicate how the stresses associated with 
vertical motion and pitching can be combined with those 
associated with rolling due to vertical gusts, particularly in 
the case of a large rigid airplane. (A small airplane, flexible 
or rigid, does not roll as a result of the action of vertical 
gusts, and for the large flexible airplane it is more convenient 
to consider rolling motion simultaneously with the other 
motions, so that the superposition is effected automatically 
in the process of obtaining the required transfer functions.) 
For all airplanes the effect of side gusts can then be taken 
into account, if .isotropy is assumed, by adding the stress 
spectra directly. 

In this section the instantaneous stress at a given part of 
the wing will be assumed to be proportional to the instan- 
taneous bending moment at that section of the wing, so that 
consideration can be confined to this bending moment; also, 
the airplane will be considered to be free to move in only 
two degrees of freedom, namely, vertical motion and rolling. 
The extension of the following arguments to stresses which 
depend on the vertical shear and the torque as well, and the 
inclusion of pitching as an additional degree of freedom, can 
be effected readily and will not alter the conclusion reached 
here. 

"When rolling motions of the airplane are taken into 
account, the bending moment due to the aerodynamic forces 
associated with rate of roll and to the inertia load associated 
with rolling acceleration can be expressed in terms of the 
rolling moment by a superposition integral. In turn, the 
rolling moment can be expressed in terms of a superposition 
integral involving the instantaneous gust intensities along 
the span, the influence function y'(y) used in the fust part 
of the report, and an associated response function h M (t ) . 
Hence, equation (52) is modified by the inclusion of a third 
term and becomes 

(y)]w(U (t—h),y)dy 

H the spectrum for M B is now calculated in the manner 
used previously, the following result is obtained: 

<Pm b (u) = \S L (o}) |V^ t (w) + |flji(w) | V»,(«) + 

2R{Hi(ci))Hjr*(o>) } <p Wf ^ (w) + |Hh(<xi) i Vr <3 (*>) 4- 

}ip a ^(o>) ( 69 ) 

where is the Fourier transform of h' M {t), and where 
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<Pw C3 (u) and <Pw Ci (u) can be obtained from equation (12) with 
the following values of r 3 and 1+ respectively: 

r a (v)^J m y'(y)y’(3/+v)dy 

r<(>?)=jJ_ w 2 ‘ y'(y)y^ a (3/+v)dy 

(Tlio function r 3 (i 7 ) is four times the function r'(ij) con- 
sidered previously in connection with the rolling moment. 
This factor of four must bo taken into account in h'n(t) .) 

The spectrum ip Mg may be considered to consist of two 
parts: The first three terms listed in equation (59), but 
excluding the contribution of y Ma to <p af (see eq. (55)), 

represent the contribution of the sy mm etric parts of 
the instantaneous gust distributions or the contribution 
associated with symmetric motion; the last two terms 
and the contribution of y Ma to <p„ represent the contribution 

of the antisymmetric parts of the instantaneous gust distri- 
butions or the contribution associated with rolling motion. 
Therefore, the power spectrum of the stress due to gusts 
and combined symmetric and rolling motion resulting from 
the gusts is the sum of the two power spectra (that for the 
gusts and symmetric motion alone and that for the gusts 
and rolling motion alone), provided the direct contribution 
of the gusts is split up into a symmetric and antisymmetric 
part and each is taken into account only once, in connection 
with the appropriate type of motion. The cross-correlation 
terms between the symmetric and antisymmetric contri- 
butions to the stress can be shown to involve integrals which 
contain products of symmetric and antisymmetric influence 
functions and, hence, vanish, so that the cross correlations 
are zero ; hence, the two parts of the combined spectrum ore 
statistically independent and, therefore, directly additive. 

DYNAMICS OF THE. SMALL FLEXIBLE AIRPLANE 

The purpose of this part of the report is to consider the 
transfer functions relating the stresses at various points of a 
small flexible airplane to the vertical gusts which cause them. 
The longitudinal and lateral degrees of freedom are still sepa- 
rable, and only the longitudinal degrees will be considered; 
the lateral degrees, which are involved in analysis of the re- 
sponse to side gusts acting on the vertical tail, can be analyzed 
in the same way. Therefore, for this case only one-half of 
the wing need be considered as a result of the symmetry (or 
antisymmetry, in the case of the lateral degrees of freedom) 
of the problem. 

The method which is outlined in this part consists in an 
extension of the numerical-integration method of aeroelastic 
analysis described in reference 14 to sinusoidal motions of the 
airplane. This extension takes into account the facts that 
the aerodynamic forces now have out-of-phase as well as in- 
pliaso parts and that vertical deflections must now be calcu- 
lated separately because the structural deformations can no 
longer be characterized by angle-of-attack changes. Also, 


the “rigid-body” degrees of freedom (vertical and pitc hin g 
motion of the airplane as a whole and structural deformations 
of the tail) are now taken into account. 

The result is a set of linear algebraic equations (which serve 
as an approximation to the ordinary differential equations 
that characterize the problem) for the airplane motions and 
deformations in terms of the applied aerodynamic forces di- 
rectly due to gusts. The desired transfer functions can then 
be obtained from solutions of these equations at various fre- 
quencies, and the power spectra of the stresses are given by 
the product of the absolute square of these transfer functions 
and the point power spectrum of the vertical component of 
turbulence. The same equations may be used to calculate 
with little additional effort certain aeroelastic effects, such as 
the static aeroelastic deformations and the flutter speed, 
which are usually obtained in separate analyses. 

LOADS APPLIED TO THE WING 

The loads applied to the wing stem from three sources: 
The aerodynamic loads directly due to the action of the . 
gusts, the aerodynamic loads due to the motions of the air- 
plane, and the inertia loads. 

The lift and pitching moment (about the elastic axis) per 
unit span on a two-dimensional airfoil undergoing sinusoidal 
angle-of-attack changes and vertical motions in incom- 
pressible flow are (see ref. 15) 

^] + 

a ~ ik w?r 

\j 2 k ‘ r/2 

The terms multiplied by C(k) are referred to as the circulatory 
terms because they are calculated from the bound and shed 
vorticity, and the others axe referred to as the potential terms. 
The potential terms are in the nature of additional-apparent- 
mass effects, and all those that involve P are usually treated 
together with the inertia forces rather than with the aero- 
dynamic forces. For compressible flow, however, the forces 
are calculated in a different manner, and the division of the 
forces into circulatory and potential parts then has little 
meaning. Consequently, in order to facilitate the extension 
of this analysis to compressible flow, this distinction will not 
be made herein. 

The aerodynamic forces are therefore written as 
l a =2irqc ||<7i(fc)a+ci,(&) 
m a = 2*qc* [&(*)«+&(*) 



i a =2 wye? JyC(k)e i (1+2 e 4 ik) 


l a —2irgc-$^C(k) (1+2 e t ik)a—ik 
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so that for incompressible flow 


am = a+ 2 e 4 ik)cm + l 4-em 

£ 

am=-vccm + f 

am = (1+2 eMeicm -f ^-+(^+^) a 

am^-i^cm-T 1 * 

£ 

In order to calculate the lift at a given point of a wing of 
finite span an appropriate Green’s function is required. An 
approximation to this function based on a reciprocity theorem 
of linearized lifting-surface theory, is given in the first part 
of this report; as used for the computations of that part, this 
function implies the assumptions that the spanwise distri- 
bution of the lift for oscillations of the wing as a whole is 
substantially invariant with frequency and that this distri- 
bution is the same in direct as in reverse flow. Neither of 
these assumptions is essential to the analysis but both, and 
particularly the first, simplify it considerably. With these 
assumptions, the desired lift distribution is then given by 
expressions of the form of equations (32a) and (34). - 
For the present purpose, however, a set of aerodynamic 
influence coefficients is required, rather than influence 
functions. Such a set of coefficients, based on the same 
ideas, can be obtained readily by the techniques used in 
references 9 and 10. The result may be expressed as follows: 

{n^ajq-^am isiw+w wa{p} } ( 62 a) 

where the aerodynamic-influence-coefficient matrix [Q] is 
defined by 


an appropriate. Green’s function could be estimated if the 
lift distribution for wings with parabolic camber were known ; 
however, such lift distributions do not appear to have been 
calculated for wings of finite span in unsteady flow. In fact, 
relatively little is known about moment distributions even 
in steady flow. However, the available information in- 
dicates that the local center of pressure does not appear to 
be very sensitive to the lift distribution. (See ref. 9, for 
instance.) It will be assumed that this is also true in un- 
steady flow at a given frequency, and that, furthermore, 
these centers of pressure are given by two-dimensional 
theory. With this assumption the moment can be written 
as 


{m T } a =C La cg.^C 3 m[c 1 [«]{«}+ G,(Je)[c\ ^{jj ^ } (62b) 


The lift and moment distribution due to the gust can be 
calculated in a similar manner. In the following analysis, 
the magnitude of the gust intensity is considered to bo 
unity, and the longitudinal reference point is the intersection 
of the elastic axis and the wing root, so that the instantaneous 
gust intensity at any station y is 


w=e 


-Italian A 
*/2 


(63) 


With this function w, 

C La cq<t>m 


{ l }i 


U 




, , C La cq4>m , imif , 

{ m T } t = [e,c\ [Q}{w] 


(64) 


where <f>(k ) is the Sears function, as before. 


[Q] = (l— -EDM { l } HI b+f-STv] 


Finally, the inertia loads are 


where, in turn, {l} is a unit column matrix, and |/J is a row of 
integrating coefficients suitable for integrating a continuous 
function for a range of its argument from 0 to 1. Thus, 
for instance, if n equidistant points on the semispan are 
considered, and n is odd, then according to Simpson’s rule 


l i =—rhz J rme^ca 
m r< = — (r*+ (e*c) ! ) m a+ me^cz 


, r , 1 1424 41 

n~l [_3’ 3* 3’ 3* 3* 3_ 


Very little information is available concerning the span- 
wise distributions of the pitching moment on wings of finite 
span in unsteady flow. By means of the reciprocity theorems 
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/ ' >=, 


wing can thus be 

written, in summary, 

2 c La d 2 (k)[Q\+ 

c La ci\mo\~ 
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( 66 ) 


C L l<i>{k) 
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Uc][<3] 


LOADS APPLIED TO THE TAIL 


The loads applied to the tail are similar in nature to those applied to the wing, but the tail experiences additional loads as 
a result of the downwash produced by the lift on the wing. Again, little is known about the downwash in unsteady flow, and 
even in steady flow the downwash cannot be predicted accurately because of boundary-layer effects on the fuselage and the 
wing root. Consequently, even in steady-flow analyses experimental results are usually relied upon. 

In the following analysis, the assumption is therefore made that experimental results are available for steady flow, in the 
form of the downwash derivative de/da. In order to determine the attenuation of this value with frequency, the results of 
the analysis of reference 16 will bo used. These results indicate that the time variation of the tail lift due to the downwash 
caused by the wing lift which results from a unit jump in the wing angle of attack can be approximated by an immediate jump 
in the tail lift of — 0.16 of the steady-state value and another jump to the steady-state value after the time required to travel 
the distance from the 45-percent-chord point of the wing to the quarter-chord point of the tail plus another eighth of the 
chord length. Hence, for sinusoidal angle-of-attack changes the tail lift due to downwash is 


(£0<,.=-<V‘S‘ [-0.16+1.16^ K^fO] [ (1+2 e,ik)a-ik (67) 

where x,' is the distance from the intersection of the elastic axis and the wing root (assumed for this purpose to be at the 45 
percent point on tho root chord) to the aerodynamic center of the tail. As pointed out in reference 16, this approximation is 
valid only for &<0.35; however, this range is adequate for the present purpose. 

Similarly, the downwash associated with the wing lift due to sinusoidal gusts gives rise to a tail lift which, within this 
approximation, is 


(I,.),— 5. § | 


-O.WJt ~ 


C W T 


( 68 ) 


(The additional lag represents the time required to travel 
the 0.6 root semichords from the 45-percent-chord point of 
the root, which is the reference point for the gusts, to the 
76-percent-chord point of the wing root, which is assumed 
to be the point governing the lift at the wing root, inasmuch 
as it is the centroid of the influence function for the chordwise 
pressure distribution.) 

The other aerodynamic forces are those due to the motions 
of the airplane, those due to the tail deformations, and those 
directly due to the gusts. On the basis of the assumptions 
made in the preceding section, these forces are 

L t tt =0 Lat i,S t [a(i)« r +C,(« =^+C 1 (k)Aa+C i {k)^ (69) 



where 
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and a, and z r are the angle of attack and vertical displacement 
of the airplane at the wing root. 

The inertia load on the tail is 


or 


Li f = — m t (z r -\-Az—Xi'a f ) 



(Zr+AZ—x/Or) 


(71) 



168 


REPORT 1345 — NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


Here the center of gravity of the tail has been assumed to 
coincide with its aerodyn ami c center; in order to remove this 
assumption it is necessary only to add (or subtract) the 
distance between the two to x/ in equation (71). 

The normal forces on the tail can then be summarized as 
follows: 

L^qiF^Zr+FtWar+Ftik) Aa+F t (k)Az ±F B (k}w t ] (72) 
where 

W=2^|f {&(*)+ 
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and 


w t =e ^ 


The pitching moments corresponding to these normal 
forces can be obtained in the manner employed for the wing. 
However, inasmuch as the tail chord is usually small com- 
pared with the fuselage length, the travel of the center of 
pressure of the tail is small compared with the length x,'. 
Hence, the center of pressure will be assumed to remain at the 
aerodynamic center of the tail, and the pitching moments 
are then — »/ times the corresponding forces, so that 

M t =— gx,'[F 1 (k)z T +F 2 (k)a r +F z (k)Aa+F i (k)Az+F 6 (k)wi\ 

(73) 

WING AND TAIL DEFORMATIONS 

The wing deformations may be calculated either from 
structural influence coefficients or from the bending and 
torsion stiffnesses of the wing used in conjunction with 
simple beam theory. The latter approach will be followed 
here, based on the method of reference 12. • 

The bending and torsion moments on the wing structure 
may be obtained by integrating the applied loads. If numer- 
ical methods are employed to perform these integrations the 
results may be written as follows: 



■(Ml 

cos A 


HU 


-sin A - [I] 



r (^) ,[II "ihr/lj 


[ 0 ] 

Similarly, the deformations are 
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(75) 


where the integrating matrices (see ref. 12) perform the following operations: 

[H]{/}oJ[JV(rww 
mmo P P'firw'df 

Jo Jo 
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These integrating matrices may he based on the trapezoidal 
rule, Simpson’s rule, or any other numerical method; the 
intervals chosen for | (0^£^1) need not be of constant 
width unless a specific rule demands a uniform spacing. 

The structural deformations of the wing may then be 
written in terms of the applied loads as 


mobile (a r =z r =0), they would be sufficient to calculate all 
unknown quantities. However, if the fuselage is free to 
move, two additional equations are required to obtain the 
two additional unknown quantities a, and z r . These addi- 
tional equations are those expressing the dynamic equilibrium 
of the forces on the fuselage, namely, 


{ z-z r 
a—a T 



_© 

© 


{w} 


(76) 


where the submatrices [®], [<§)], [®], [®], [©], and [(f)] 
designate, respectively, the four quadrants and two halves 
of the matrices, obtained by postmultiplying the square 
matrix of equation (75) by the square and the rectangular 
matrix of equation (66) , respectively. 

For the purpose of the present analysis, which is concerned 
primarily with the wing stresses, the tail deformations are 
treated by including only the vertical displacement and 
angle-of-attack change of the tail as a whole due to the tail 
load. These quantities may be obtained from a static test 
which consists in applying a concentrated normal load at 
the aerodynamic center of the tail and measuring these 
deformations. They may also be obtained from a vibration 
test in which the deflection curve of the rear part of the 
fuselage in the lowest vertical-bonding mode is measured; 
in this case the desired spring constants can be deduced from 
the frequency relation of a simple mass oscillator in terms 
of the measured frequency and of the mass of the empennage 
(including the part of the fuselage which may be considered 
to move with the empennage) . In the absence of such tests 
these constants may bo calculated in an analogous manner. 

The tail deformations may then be written as 


so that, also, 


A a= — KiL t \ 
ts.z=K i L l J 


. K 3 . 

Az= — A a 

Hi 


(77) 

(78) 


Inasmuch as these deformations are not independent of each 
other, only one need be retained in the analysis. Therefore, 
if As is eliminated by means of equation (78), Aa can then 
be obtained in a form similar to that used for the wing 
deformations in equation (76): 

Aa = — qKi Fi (k)z r +F 2 (k)a r +^F 3 (k) — 

§ W]a«^ -gK,F 6 {k)w t (79) 

EQUATIONS OF MOTION 


L w -\-L—m f (z r — Ar 0^=0 
M w —x t L t — (r/+ (Ar) 5 ) m^+Ar m / z r =0 


(80) 


where the wing lift and pitching moment can be expressed in 
terms of the lift and moment distributions l and m as 


M„=\ -2 (£)’ tan A [IIJ, 7.\ [I] f l \ 


(81) 


or 


£*= 21 ®] 


M a =q[@ J 



-fg[©J{w} 

+2l@J{w) 




(82) 


where [@], [©J, [@J, and [@J are the rows obtained by post- 
multiplying the rows of equation (81) by the square and 
rectangular matrices of equation (66). In equation (80) the 
fuselage lift and moment have been neglected; they can 
easily be expressed in terms of a r and z r and included, if 
desired. 

The equations for the tail deformation Aa and those for 
the overall normal force and pitching moment can be com- 
bined with equation (76) as follows: For the sake of definite- 
ness it will be assumed that n stations on the wing are 
considered, including the one at the root, so that there are 
271+ 1 unknown quantities, and that in the column matrices 
defining applied loads, deformations, and so on, the values 
at the root of the wing are written at the top: 





Equations (76) and (79) are equations of motion inasmuch 
as they describe balances of aerodynamic, structural, and 
inertia forces. In fact, if the airplane fuselage were im- 


The first and ( 71 + 1 )th equations of the system defined by 
equation (76) express only the trivial fact that the structural 
deformation at the wing root is zero. They are replaced by 
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equations (82) and (79), which are adjoined to the system, to yield the combined equation of motion 
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(83a) 


If the square matrix on the right side of this equation is 
designated by [4], the rectangular matrix by [B\, and the 
quasi-unit matrix on the left side of the equation by [1'], the 
equation can also be written as 







— F s (k) 




(83b) 


The matrices [©] to [©] are the same as the matrices 
[®] to [©] of equation (76), except that the first rows of 
the latter, which are all zero, are replaced as follows: 


First row 
of— 

Replacement 

Quantity added to the leading element 

© 

First half of l@J 

2\2/ 

© 

Second half of l(u)J 

«*> ,/T> 

2\2 ) 

© 

First half of l@J 

-xJfiWr-ZgkV 

1(1) 

© 

Second half of l@J 

xF,(k)+ m ' [r ''+^ ] V 

1(1) 

© 

l@J 

0 

© 

[®] 

0 


Also, the elements of the last rows and columns of the mat- 
rices [A] and [15] are zero, except for the following: 

A x . an+i =F,(k)-~?F 4 (k) 

A n+ u as+1 = -x, [>,(£) -fpw] 

■rijn+i, i= K\F i{k) 


■4sii+i, *+i K\F a (k) 


Bn+l.n+l — X,Fi(k) 

B2X+1. *+i = — K\F r,{k) 

SOLUTION OF THE EQUATIONS OF MOTION 

For the purpose of calculating the desired transfer func- 
tions, equation (83b) may be solved directly for a given 
value of q as a set of linear algebraic equations with coeffi- 
cients given by the matrix [l']—ff[.A] and with “knowns” 

given by the column matrix q[B] {^,j K ( Where antl 

™ are functions of kj. The result is a column matrix of 

the unknown amplitudes of the motions of the airplane. If 
this column is calculated for several values of k in the range 
of interest, these amplitudes, considered as functions of k, 
are transfer functions from the gust to the motions. 

This column matrix can be substituted into equation (66) 

and the resulting column matrix {JJ substituted into equa- 
tion (74), to yield a column matrix of bending and twisting 
moments which again, considered as a function of k, repre- 
sents transfer functions from the gust to these moments. 
A set of transfer functions for the vertical shear could bo 
calculated similarlv from the relation 




The stress at any point of the structure can be assumed 
to be given by a linear superposition of the bending moment, 
twisting moment, and vertical shear at the given station, if 
elementary beam theoty is used. If elementary beam theory 
cannot be used because of the interaction of bending and 
torsion stresses or because of shear lag, the stress at a given 
point can be expressed as a linear superposition of moments 
and shears at other stations as well as the given station. In 
either case, the transfer function for the given stress is then 
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produced by the same linear superposition of the transfer 
functions for the corresponding moments and shears. 

It may be noted that, at zero frequency, solution of equa- 
tion (83b) yields the static aeroelastic deformations and thus 
permits the calculation of the changes in the lift distribution 
and the shift of the aerodynamic center that result from 
static aeroelastic action. Also, inasmuch as this equation 
completely describes the dynamic behavior of the airplane, 
the speeds at which aeroelastic instability phenomena occur 
can be calculated from it, although such calculations are 
beyond the scope of this report. Suffice it to point out that 
for such a calculation the degrees of freedom of the body 
must bo eliminated first, as a result of the way in which 
the problem has been set up. This elimination can be 
effected readily by considering the first and (n-(-l)th rows 
of [A], but with A u , A, *+i, A nHwU and A n+hn +i replaced by 
0, If these rows are premultiplied by 

pin "I -1 

L-^x+l.l -d.n+ 1 , h +1 J 


and used as the first and («-)-l)th rows of a matrix which 
is otherwise a unit matrix, and if this resulting matrix is 
referred to as [1"], then equation (83b) can be written for 
this homogeneous case as 


m in~g.[A][i"]] 



= { 0 } 


The products [1'] [1"] and [ri] [1"] will now have two null 

rows and columns each, which correspond to Zr and « r . If 

these rows and columns are deleted and z, and a r are deleted 

/ 



, yielding a column 



, the remaining mat- 


rices are nonsingular, so that they can be inverted and the 
preceding equation can be mitten as 


[[.E]- 2 [D(A:)]] 


{£}'- 


(84) 


where [E] is the identity matrix, and 

[7?(A:)] ^ [ [PfU 77 !] _1 CU] [I'T ] 


The horizontal braces designate the fact that the null rows 
and columns have been deleted. 

Equation (84) is in the canonical form for the calculation 
of eigenvalues. If k is set equal to zero and the eigenvalues 
of [0(0)] 11X0 calculated by iteration, expansion of the de- 
terminant, or any other suitable method, the lowest real 
and positive one represents the value of the dynamic pres- 
sure at divergence. For swept wings the value lowest in 
absolute magnitude is usually negative and is therefore of 
no practical significance, although it is often used as an 
index of the aeroelastic behavior of the airplane. 

This calculation can be repeated for various positive values 


of Jc, the first few eigenvalues being obtained for each. The 
results, which will generally be complex, can be plotted 
against k. When any of the eigenvalues becomes purely 
real, it represents a dynamic pressure at flutter, and the 
corresponding value of k represents the reduced frequency 
at flutter. (This statement is true only if the structural 
damping is zero; such damping effects can easily be included, 
but the details of the process are beyond the scope of this 
report.) 

DYNAMICS OF THE LARGE FLEXIBLE AIRPLANE 

For the large flexible airplane the fundamental proposition 
of power spectral analysis, that the output power spectrum 
of a system is the product of the absolute square of the 
transfer function and the input power spectrum, is no longer 
valid if the input is considered to be the gust intensity at a 
point. Nor can the output power spectrum be expressed 
directly in terms of an effective input spectrum, as in the 
case of the rigid airplane, where this simplification resulted 
from the assumption that the indicial-response function was 
expressible as the product of a function of time alone and a 
function of distance along the span alone. 

In the first section of this part of the report the statistical 
problems involved in an analysis of the response of a large 
flexible airplane are considered. The nature of the general- 
ized transfer functions required for this purpose is described, 
and the means whereby they are combined with the input 
spectrum are indicated. 

The second section is concerned with an extension of the 
method outlined in the preceding part to the case of the large 
airplane. Although fundamentally the dynamic aspects 
of the problem are unchanged, and although the longitudinal 
and lateral degrees of freedom can still be separated, a direct 
application of the approaches outlined in the preceding parts 
of the report to the large airplane requires consideration of 
the entire wing, rather than only one half of the wing. 
Little additional computing time is then required to treat 
the lateral and longitudinal degrees of freedom simul- 
taneously, and the necessity of combining the results of two 
separate analyses is obviated. However, attention can 
still be confined to one half of the wing by using the technique 
outlined in the discussion following equation (55) and in the 
section headed “Combination of the Results Obtained From 
Analyses of the Longitudinal and Lateral Degrees of Free- 
dom.” Basically this technique consists in splitting the 
influence functions of concern into symmetric and anti- 
symmetric parts and using one part for an analysis involving 
the longitudinal degrees, and the other in an analysis in- 
volving the lateral degrees. If this approach is adopted, 
separate analysis of longitudinal and lateral degrees of 
freedom is still preferable. 

EXTENSION OF THE STATISTICAL APPROACH 

The power spectrum and, hence, the mean-square values 
of the responses such 'as the stress at a given point on the 
wing of a large flexible airplane due to flight through turbu- 
lent air can be calculated in several ways. Perhaps the most 
direct of these consists in using the basic approach outlined 
in the first part of this paper and starting with an expression 
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for the instantaneous value of the stress in terms of a suitable 
indicial-response influence function, namely, 

cm 

cr(t)= K{U,y)w(lJ(t—U) ,y) dy dU (85) 

J —co J —6/2 

The function K° c (t,y) is, as before, an indicial-response 
influence function or Green’s function for the partial differen- 
tial equation. It relates the stress as a function of space 
and time coordinates to the applied loads, which are also 
functions of space and time coordinates. The essential 


difference between the problem considered in this section 
and those treated in the first two parts of the report is that 
this function can no longer be expressed as a product of a 
function of t alone and y alone. In some cases — in a modal 
approach, for instance — it may be expressible as a sum of 
several such functions, and then the approach used previously 
for the bending moment of a large rigid airplane free to move 
vertically may be adopted, as has been done in reference 13. 
However, in this section the case is considered in which 
even this simplification cannot be made. 


The correlation function for <j can be calculated directly from equation (85), and for the case of axisymmetric turbulence 


it is 


Mt)= f f f f K{ti,y^h%(t it y^ {r}U i {T-\-ti—t^ i +(y 1 —y^) dy x dy 2 dt x dt a 

%) — ® ti — <x> \) — 6/2 %) “6/2 


( 86 ) 


Hence, the power spectrum of cr can be obtained by calculating the Fourier transform of this function and is 

J 'bl 2 C6/2 

H~\u,y l )Hz(w,y£vu(<>>,\yi—yi\)dy 1 dy i (87) 

-6/2 J -6/2 

where ^(w,*?) is the two-dimensional spectrum considered previously, H^(co,y) is the Fourier transform of h*(i,y) with respect 
to time, and the asterisk designates a complex conjugate, as before. Thus, the function H^{w,y) is a transfor function from 
sinusoidal vertical gusts (of width dy) impinging on the wing at a given station y on the wing to the stress a-, or a Green’s func- 
tion for the ordinary differential equation (with the quantity u as a parameter) that relates the stress amplitude os a function 
of the space coordinate y to the amplitude of the applied sinusoidal gusts. 

The term dTf(ai,y{) IT c (w,yO in equation (87) is complex; however, the imaginary part can be ignored because it contributes 
nothing to the integral as a result of the fact that <p a depends only on |y 2 — l/i|. 

In a maimer analogous to the one employed in the first part of this report, the double integral in equation (87) can bo 
evaluated by introducing the function which takes the place of the function T(t?) used previously and is defined by an 

autoconvolution of H^(«,i?) : 


( 88 ) 


~ P6/2 

S’( 03 ,t {)— 2 K. { (co, y) H*(u ! y+rj)}dy 

J ~bfi 

where R{ } designates the real part. Hence, 

X b ~ 

Hcfarfipufa, 77) dr] 

Another approach consists in using the double Fourier transform $„(X) of the axisymmetric input correlation defined in 
the first part of this report. In terms of this function the correlation function &(r) can be written as 

1 ✓»<» nm nbtl 

Mr) =7 K {h,ydK (ti,yj (Vx^+Xj*) d\ x d\dy\dy,flhdU 

— coi/ — an/ —6/2*7 — 6/2*/ — a>*7 — co 

=4 r r e ^|^ ( x 1 ,x 2 )| 2 $tt (vv+x 2 *) dx^x, 

— CO x) — co [ ! 


where the function 

* c° cm 

I?“(Xi,Xj)= h%(t,y)e- i< * ut+ w dydt 

*/ — 00 */ — 6/2 

represents the Fourier transform (with respect to y) of the function H^(o>,y), that is, 


^(X 1 ,X a )= f" dy 

J-b/S 


( 80 ) 
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Hence, 




$ u (VV+V) d\idki dr 


-2f£Kp' x )I , *'(V© ,+x, )‘ !x 


(90) 


With the approach outlined in the next section, the function E%(o>,y) can be calculated either directly or indirectly, by 
first calculating the transfer function from the gust to the lift distribution and then the transfer function from the lift distri- 
bution to the stress. For the indirect method, 


J i/l 

E l a (.<->, y) E ” (w) G(v,y)dy 

-V a 


where the function E^(cS) G(t],y) is the influence function for the local lift or the Green’s function for the aerodynamic 
problem involving sinusoidal gusts considered in the first part of this report; the symbols y and y in G(y,y) are inter- 
changed, however, so that the function now defines the contribution of a gust at station y to the lift at station y. 
The transfer function Hl(a),y) relates the (sinusoidal) stress at the given point to a unit concentrated (sinusoidal) normal 
force acting at station y. 

With this indirect method, the power spectrum for the stress at a given point can be calculated by starting with the 
power spectrum for the lift distribution calculated in the first part of the report. For this approach <r(t) may he 
writton as 


so that 

and 

where 


/•<= /»S /2 

<r(t)= K{ti,y)l{t—ti,y)dydti 

«/ — »*/ — 6/2 

J co r»co *t ,/ 2 pfc/j 

- K{ti,ydK{ti,yi)'l'i(j+tx—ti,yi,yddy l dydtidt i 

- co */ — co — 6 / 2*7 — 6/2 

J &/2 rt >/2 

IP’ (co,y,) E‘ (u,y a ) <p, (cu,j/i ,yjdyi dy? 

-vs J-m 


(91) 


The Fourier transfoim of this correlation function is, then, 

<pi (u,yi,y 2 ) = |H’“(co) | V» e (<■>) 

(See the section headed “Generalized Aerodynamic Influence Functions in Unsteady Flow.”) Hence, 

nn s nt>r a 

^(co) = |ZZf(co) I 2 fl 1 ; (co ; y i) (a>,y 2 ) <p v (to) dy x dy 2 

J-mJ-m 

If <Pw t (a) is given by equation (35), the double integral can be expressed in terms of single integrals as follows: 

^(w)=|j??(<o)| 2 £(1— K)*\ G(<d) | : i <p at (b>)+2K(\—K)'R { G(cS) G'*(o}) <?"(&>, i?)¥>»(w,i?)di?J 


(92) 


(93) 


( 94 ) 


where 


Jin 

Ei(a>,y)y(y)dy 

-W 2 

Ji/i 

HIM**- (o>,yh(y)dy 
-m 
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as before, R{ } designates the real part, andl{ } designates the 
imaginary part. 

Equations (88), (90), and (93) thus represent three 
methods of obtaining the power spectrum of the given stress. 
One requires a transfer function from the local gust intensities 
directly to the stress, an autoconvolution 'of this transfer 
function, and the spectrum of turbulence defined by equation 
(16) ; the second requires a two-dimensional spectrum of the 
turbulence defined by equations (17a) and (17b) and a 
Fourier transform of the aforementioned transfer function 
with respect to y; the third utilizes an autoconvolution of 
the transfer functions from local concentrated loads to the 
stress and the spectrum for local lifts due to vertical gusts 
calculated in the first part of the report. The choice of 
approach depends to some extent on the information avail- 
able, but is largely a matter of individual preference. 

The analysis in this section has been based on the premise 
that both halves of the wing would be treated simultaneously. 
Attention can be confined to one wing by using only the 
symmetric parts of the influence functions in an analysis 
involving the longitudinal degrees of freedom, and only the 
antisymmetric parts in an analysis involving the lateral 
degrees of freedom. The symmetric part of an influence 
function for a unit concentrated load or gust acting at sta- 
tion y is the response function for two loads or gusts of % unit 
intensity acting at stations y and —y, respectively. Simi- 
larly, the antisymmetric part is the response for a load or 
gust of K unit intensity acting at station y and an equal and 
opposite load or gust acting at station —y. If this approach 
is chosen, some of the integral expressions given in this 
section assume slightly simpler forms. For instance, equa- 
tion (89) can then be -written for the symmetric part as 

* rt>n 

HJ(X 1 ,^ a )=2j cos (Kill ,y)dy 

and for the antisymmetric part as 

a rm 

(X! ,X 2 ) = — 2i sin \ 3 yH?(\ 1 U,y)dy 

Jo 

CALCULATION OF THE REQUIRED TRANSFER FUNCTIONS 

Depending on which of the methods outlined in the pre- 
ceding section is used, one of two types of transfer functions 
is required — either the one from local gusts to the stress of 
interest, or the one from the local lift to that stress. Both 
of these functions differ in several respects from those con- 
sidered in the preceding part of this report. 

For either type of transfer function the structural de- 
formations of the wing under concentrated loads as well as 
under distributed loads will be required, so that the numeri- 
cal-integration schemes used in the preceding part have to be 
modified to some extent. Also, it may now be preferable to 
consider the entire wing (rather than the one semispan con- 
sidered in the case of the small flexible airplane), so that the 
various transfer functions are asymmetric. "With the degrees 
of freedom considered in the analysis thus doubled, it becomes 


preferable to include one additional degree of freedom, roll, 
rather than to perform two separate analyses for symmetric 
and antisymmetric motions with, respectively, one and two 
less degrees of freedom, and then to combine the results. 
However, if the alternative approach of splitting up the 
influence functions into symmetric and antisymmetric parts 
for use in two separate analyses is adopted, only one semispan 
need be considered, and the results presented in this section 
can then be simplified to a large extent. 

Before discussing the modifications required to extend the 
dynamic analysis outlined in the preceding part of the report 
to the large airplane, it might be pointed out that chordwiso 
deformations (deformations parallel to the chord) will again 
be ignored. Again, they can readily be included by a 
straightforward extension of the approach used here if it is 
felt that they may be significant in any given case. If they 
are included, however, yawing and possibly also sideslipping 
motions can probably no longer be ignored, because they 
may give rise to large forces in the chordwise direction. If 
the entire wing is treated, these two additional degrees of 
freedom can readily be included, and all longitudinal and 
lateral degrees of freedom are then treated simultaneously; 
if two separate analyses are performed for the longitudinal 
and lateral degrees of freedom, the symmetric and anti- 
symmetric parts of the influence functions being used and 
only one semispan being treated, these additional degrees of 
freedom enter only into the lateral analysis. 

The structural deformation due to local (concentrated) 
loads can be obtained in several ways. If measured influence 
coefficients are used, they pertain precisely to such loads and, 
in fact, must be modified before they can be used for dis- 
tributed loads (see ref. 12) so that it is necessary only to use 
the unmodified coefficients. 

If the deformations are to be calculated in a manner 
similar to that employed in the preceding part, the inte- 
grating matrices must be replaced as follows: 


where 


and 


where 



W=1 
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The factor of % for l vv ' constitutes an approximation which 
implies fairing through a discontinuity. If this approxi- 
mation is to be avoided, the deflections due to unit concen- 
trated loads (the structural influence coefficients) can be 
calculated directly from simple beam theory, in -which 
case the limits of integration take care of the discontinu- 
ities. Thus, for instance, for an unswept wing, the normal 
deflection and twist at y v due to a unit concentrated load 
and torque, respectively, at y s are 


Jo ’J EI{y) dydy ' (2/p - Vq) 

2m= /o Jo w§) dvdy,Jr(:!Jv ~ y,) i Q w§) dy 



(Vr ^ Vt) 


(.Vc < Vv) 


The concentrated loads under consideration arise as 
follows: For the transfer functions relating local lifts to 
the desired stress, the local lifts may be considered to be 
concentrated loads of unit magnitude, associated with 
concentrated torques of magnitude e,c. Equation (83b) 
can then be -written as (see also eq. (76)) 

[ [i ']- S M][£]=m[i^ 5 J <w>. 


where (to] is a diagonal matrix of the values of w defined by 
equation (63), and where the matrix [22'] represents either 
the four influence-coefficient matrices for z and a due to 
concentrated loads and torques, or the square matrix of 
equation (76) with modified integrating matrices, as dis- 
cussed in the preceding paragraphs. 

It may be noted that equation (95) now represents not 
one set of simultaneous equations but several, all having 
the same coefficients but different sets of knowns (as defined 
by the columns of the matrix on the right side) and, 

hence, different sets of unknowns ^the col umns of the 


matrix 



This situation is due to the fact that the 


functions under consideration are, in effect, the responses of 
the airplane as a whole to sinusoidally varying concentrated 
loads and are different for each location of the applied load. 

Once this equation has been modified to take into account 
the overall body motions and tail deflections (as explained 
in the preceding part) as well as the rolling motions (as 
explained in the following paragraphs), it can be solved to 


yield the unknown values of 



From these values the 


bonding and twisting moments, as well as the vertical shears, 
can be calculated and added to those due to the concentrated 


loads. Whan combined linearly as required for the desired 
stress, these moments and shears yield the desired transfer 
functions Hl(ca,y). 

If the transfer function directly from the local gusts to 
the desired stress is to be determined, the response of the 
airplane to the lift distribution induced by a sinusoidal gust 
of width dy acting at station y must be calculated. This 
lift distribution is the Green’s function considered previously. 
If it is represented by the relation 27? (o>) G(y,i]), with 6(y,tj) 
defined by the approximation given in equation (34), the 
concentrated loads arise from the delta function in that ex- 
pression. The right side of equation (95) becomes, in that case, 


a(i-K)Et 



[ y (y) y fa)] 

[ercTMyjyiv)] 



qbKHf(u) [R'] 


l y(y)l M 
.kcf[r (i/jiM- 


where [22] represents the square matrix of equation (75), 
and [R'] the one discussed in connection with equation (95). 
Again, several sets of simultaneous equations are implied. 
Their solution (after modification for overall motions and 
tail deflections) yields values of z and a from which the 
transfer functions Iff(w,y) can be calculated. 

The extension of the method of the preceding part to the 
calculation of the deformations on both wings is straight- 
forward. Essentially, distributed lifts and torques now have 
to be calculated for both wings and integrated with matrices 
which can be assembled from those used for one wing alone. 
No new problems arise in this process, so that it need not 
be discussed further. 

The inclusion of rolling motion, however, is not so straight- 
forward. One method consists in replacing all values of z 
in equation (83a) or its equivalent by z-\-9y and then 
reducing the columns involving this quantity by the following 
relation (which assumes that the new unkn own quantity, 
the roll angle 0, is listed at the end of the column): 
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An additional equation must then be joined to the set, 
namely, the equation of equilibrium in roll 

J ' 6/2 

l(y)ydy—I x e—M Dt e=0 (96) 

— 6/2 

where I x is the inertia in roll of the fuselage and empennage 
alone inasmuch as the inertia effects of the wing are included 
in l(y), and M Dt is the coefficient of damping in roll for the 
empennage. For most cases both of these contributions are 
negligible. If they are ignored, equation (96) can be written 
in matrix notation as ' 

inj{0=o 

where [IIJ is now a matrix which serves to perform the inte- 
gration required in equation (96). This condition can then 
be adjoined to the other equations of the set in the same 
manner as equations (82) were adjoined to the set in the 
preceding part. 

The result, again, is a set of simultaneous equations for z, 
a, A a, and 9 from the solution of which the desired transfer 
functions can be obtained as outlined in the preceding para- 
graphs. Also, as before, once the unknowns z T , a r , and 0 are 
eliminated from the set, the divergence and flutter speeds 
can be calculated by conventional matrix operations; these 
speeds will then pertain to an airplane free to move verti- 
cally as well as in pitch and roll and, hence, will include 
divergence and flutter speeds in antisymmetric as well as 
symmetric modes. 

DISCUSSION 

SOME IMPLICATIONS OF THE ASSUMPTIONS CONCERNING THE NATURE 
OF ATMOSPHERIC TURBULENCE 

The turbulence was assumed to be homogeneous in order 
to make the problem stationary in the statistical sense and 
thus permit the use of the mathematical techniques developed 
for such problems. In a practical sense, turbulence can be 
homogeneous only in a limited body of air. The assump- 
tion thus implies that the dimension of this body of air 
along the flight path is large compared with the distance 
traversed in the reaction time of the airplane, which in the 
case of load studies is of the order of the time to damp to 
one-half amplitude, but in the case of motion studies may 
be much larger. Obviously, the greater the body of air, 
the greater the reliability with which the loads and motions 
can be predicted (in a statistical sense) for one run through it. 
In general, turbulence at very low altitudes, which may be 
influenced significantly by the configuration of the ground, 
and turbulence in thunderstorms may not be sufficiently 
homogeneous for the purpose of this type of analysis, but 
other types of turbulence are likely to be substantially 
homogeneous over sufficiently large distances. 

Isotropy was assumed in order to permit the required two- 
dimensional correlation functions to be expressed simply in 
terms of the one-dimensional correlation functions. For 
sufficiently short wave lengths all turbulence is isotropic, but 
for long wave lengths it can be isotropic only if it is homo- 
geneous (both in the plane of the flight path and perpendicu- 
lar to it). The condition of axisymmetry, which is sufficient 


for most of the results presented herein, is less restrictive than 
isotropy inasmuch as it does not specify the variation of the 
characteristics of the turbulence in the vertical direction. 
In practical problems, if the turbulence may be assumed to 
be homogeneous, the conditions of axisymmetry and iso- 
tropy are likely to be satisfied to a sufficient extent to permit 
the use of the approach presented herein for all but very 
long wave lengths. The wave length at which it ceases to 
be valid depends on the size of the body of air under con- 
sideration, being larger for a large body. 

Taylor’s hypothesis (to the effect that a space displace- 
ment Ax along the flight path may be identified with a time 
displacement r=Az/Z7 in the gust correlation functions) im- 
plies that the variation in gust intensity that prevails along 
the flight path at any instant will remain substantially the 
same until the airplane has traversed the given body of air. 
The required correlation functions for atmospheric turbulence 
are thus in the nature of space correlation functions (rather 
than time correlation functions) and have been considered 
as such. The statistical characteristics of the turbulence 
are then independent of the speed at which it is traversed. 
Clearly, whether or not this hypothesis is valid depends on 
the flying speed of the airplane. On the basis of present 
knowledge no definite lower limiting speed can be quoted. 
However, indications are that the hypothesis is valid for 
flying speeds greater than about 100 or 200 feet per second. 
The effect of finite flying speed on the gust correlation func- 
tion can be expected to be most pronounced for large dis- 
tances, where the correlation is weak, so that the effect on 
the various spectra is likely to be small and to occur at the 
longest wave lengths, where, as previously mentioned, the 
spectrum is somewhat uncertain for other reasons as well; 
this effect is thus more likely to be significant for large than 
for small airplanes. 

The particular correlation function used herein for the cal- 
culations of the “averaged” correlation functions and spectra 
has certain theoretical shortcomings — primarily that the 
associated spectrum does not decrease rapidly enough for 
very short wave lengths. However, it does appear to be 
adequate to represent the available information concerning 
the spectra of atmospheric turbulence (see ref. 17, for in- 
stance) because the behavior at very short wave lengths is 
relatively unimportant, inasmuch as airplanes cannot re- 
spond to them, and the behavior at very long wave lengths 
is usually in doubt by reason of the nonhomogeneity of actual 
turbulence. In the intermediate range of wave lengths, this ' 
correlation function appears to be quite satisfactory. 

The parameter L* (the integral scale of turbulence) used 
herein is, for practical purposes, a largely fictitious quantity, 
inasmuch as it is proportional to the values of the gust 
spectrum for infinite wave lengths, which, in view of the 
uncertainties in the values of the spectra at large wave 
lengths, have little physical significance. Therefore, at 
present, insufficient information is available to give a value 
for L* to be used in connection with the numerical results 
calculated herein, although a value of 1,000 to 2,000 feet 
appears to be appropriate. As more information concern- 
ing the spectrum of atmospheric turbulence becomes avail- 
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able, more definite values can be deduced by fitting an 
analytical expression of the type used here to measured 
results in the range of frequencies of primary interest, and 
then using this expression as a means of obta inin g a value 
of L* by extrapolation of the measured results to infinite 
wave lengths (zero frequency). 

CONSIDERATIONS PERTINENT TO THE APPLICATION OP STATIONARY- 
RANDOM-PROCESS TECHNIQUES TO GUST-LOAD PROBLEMS 

The purpose'of this section is to point out how, in principle, 
the results of analyses of the type outlined herein may be 
used in overall load analysis and, hence, in the design of an 
airplane. 

Consideration is confined in this paper primarily to the 
power spectra of the motions and stresses of interest. ■ As 
pointed out in references 1, 2, 3, and 11, for instance, a great 
deal of statistical information of direct interest can be 
obtained from the power spectrum. For instance, if the 
random process of concern (say, the given stress as a func- 
tion of time) has a Gaussian probability distribution, the 
expected number of peaks at or beyond a given level in a 
given period of time can be calculated very simply from the 
integral of the spectrum and its second and fourth moments. 

The results obtained in this manner pertain to continued 
flight in a given body of turbulent air. They have to be 
generalized by determining the likelihood of flying through 
turbulence of the given characteristics. (See ref. 18.) The 
probability of exceeding a given stress level during the ex- 
pected life of the airplane while flying through atmospheric 
turbulence can thus be calculated in straightforward fashion. 
To this probability must then be added the probability of 
exceeding this level in maneuvers, landings, and, possibly, 
also in turbulence due to thunderstorms, because in view of 
the possibly nonhomogeneous character of turbulence in 
thunderstorms and the possibly nonlinear nature of the 
aerodynamic forces incurred while flying through them, the 
techniques used herein may not be applicable to flight 
through thunderstorms, and a separate analysis may have 
to be performed. 

Although the available information concerning atmos- 
pheric turbulence is inadequate to permit of any definite 
conclusion, the results obtainable with the approach outlined 
herein may turn out to be most significant for the prediction 
of the low- and medium-amplitude stress cycles which are 
important to fatigue studies; their validity for or contribu- 
tion to the prediction of very severe loads remains to be seen. 

CONCLUDING REMARKS 

The statistical approach to the problem of calculating the 
dynamic responses and the stresses of an airplane subjected 
to continuous random atmospheric turbulence has been 
extended in several respects; basically, only the assumptions 
of linearity, that is, of small motions and deformations, as 
well as homogeneity and axisymmetry of the turbulence are 
retained. 

The first problem considered was the effect of spanwise 
variations of the instantaneous turbulent velocities on the 
lift and moments due to turbulence. The mean-square lift 
has been shown to be reduced considerably if the span of the 


airplane is relatively large compared with the integral scale 
of turbulence. The shape of the spectrum of this lift is 
affected relatively little by spanwise variations of gust 
intensity, except at very high frequencies, if the decrease in 
the effective mean-square intensity is taken into account. 
The effect of sweep on the mean-square lift and its spectrum 
has been shown to be small for wings with a given distance 
from root to tip. 

If the variation of the instantaneous velocities is taken 
into account, the rolling moment to which the airplane is 
subjected can be calculated. The mean-square rolling 
moment has been shown to be proportional to the ratio of 
the wing span to the integral scale of turbulence for small 
values of that ratio. Similarly, expressions for the mean- 
square values and the power spectra of the local lift, the 
bending moments, and the pitching moment have been given. 
For some of these forces the required aerodynamic informa- 
tion cannot be calculated by existing methods. Therefore, 
certain approximations, based on experience with steady 
aerodynamic forces and available knowledge concerning 
unsteady forces, had to be made for the aerodynamic influ- 
ence functions in unsteady flow. 

The next problem considered was the dynamic response of 
a rigid airplane to random turbulence. This problem had 
previously been treated for the case of an airplane free to 
move only in the vertical direction and small enough so that 
variation of the turbulent velocities along the span could be 
neglected. In the present report the response of an airplane 
in three longitudinal degrees of freedom was considered; 
calculations were made which suggest that the inclusion of 
deviations from the mean horizontal motion is superfluous 
in gust-load calculations. For the remaining two longi- 
tudinal degrees of freedom, the mean-square normal and 
angular acceleration have been shown to be functions of 
only two parameters other than the mass ratio and scale 
parameter of the single-degree-of-freedom case, namely, 
dimensionless forms of the short-period frequency and of 
the time to damp to one-half amplitude. An indication is 
given of the manner in which the results obtained in connec- 
tion with the first problem can be used to extend this dy- 
namic analysis to the case in which variations of the turbulent 
velocity along the span have to be taken into account. 

The last problem treated was the dynamic response of a 
flexible airplane, including vertical motion, pitch, and, when 
necessary (as when spanwise variations in gust intensity are 
taken into account), roll. Horizontal and lateral (yawing 
and sideslipping) motions were disregarded because they do 
not generally affect the wing stresses due to vertical gusts. 
A method which represents an extension to the dynamic case 
of a numerical-integration approach to the static aeroelastic 
problem has been outlined for the analysis of the problem at 
hand. The modifications required in the basic statistical 
approach and in this method of dynamic analysis in order 
to treat the case in which spanwise variations of the gust 
intensity are important have been discussed. 

Although most of this analysis has been confined to the 
vertical component of turbulence, it has been shown that the 
simultaneous action of longitudinal, vertical, and lateral 
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gusts on the wing stresses (with due allowance for the fact 
that vertical gusts affect both the longitudinal and the lateral 
motions of the airplane) can be taken into account by simply 
adding the power spectra of the various contributions, pro- 
vided the turbulence is isotropic; the cross correlations or 
spectra have been shown to vanish either by the symmetry 
or antisymmetry of the influence functions involved or as a 
result of the statistical independence of mutually perpendicu- 
lar velocity components. 


The approach presented herein thus furnishes a founda- 
tion for the prediction of the statistical properties of the 
stress experience of a given airplane once the appropriate 
statistical characteristics of the atmosphere have been de- 
termined. 

Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Ya., November 5, 1956. 
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